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ABSTRACT 

We present a theoretical formalism for determining the structure of molecular 
clouds and the precollapse conditions in star-forming regions. The model consists of a 
pressure-bounded, self-gravitating sphere of an ideal gas that is supported by several 
distinct pressures. Since each pressure component is assumed to obey a polytropic law 
Pi(r) cx /9 7pi , we refer to these models as multi-pressure polytropes. We treat the case 
without rotation. The time evolution of one of these polytropes depends additionally 
on the adiabatic index ji of each component, which is modified to account for the 
effects of any thermal coupling to the environment of the cloud. We derive structure 
equations, as well as perturbation equations for performing a linear stability analysis. 
Special attention is given to properly representing the significant pressure components 
in molecular clouds: thermal motions, static magnetic fields, and turbulence. The 
fundamental approximation in our treatment is that the effects of turbulent motions 
in supporting a cloud against gravity can be approximated by a polytropic pressure 
component. In particular, we approximate the turbulent motions as a superposition of 
Alfven waves. We generalize the standard treatment of the stability of polytropes to 
allow for the flow of entropy in response to a perturbation, as expected for the entropy 
associated with wave pressure. In contrast to the pressure components within stars, 
the pressure components within interstellar clouds are "soft", with polytropic indexes 
7pi < 4/3 and (except for Alfven waves) adiabatic indexes j{ < 4/3. This paper focuses 
on the characteristics of adiabatic polytropes with a single pressure component that are 
near the brink of gravitational instability as a function of 7 P j and ji for 7 p j < 4/3. The 
properties of such polytropes are generally governed by the conditions at the surface. 
We obtain upper limits for the mass and size of polytropes in terms of the density and 
sound speed at the surface. The mean-to-surface density and pressure drops are limited 
to less than a factor 4 for 7 P < 1, regardless of the value of 7. The central-to-surface 
density and pressure drops in isentropic clouds (7^ = 7 p j) are also limited, but they can 
become quite large (as observed) in non-isentropic clouds, which have ji > j p i. We find 
that the motions associated with Alfven waves are somewhat less effective in supporting 
clouds than are the kinetic motions in an isothermal gas. 
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1. INTRODUCTION 

Star formation plays a crucial role in galactic evolution and in determining the structure 
of the interstellar medium. In general, star-forming regions appear to be in rough hydrostatic 
equilibrium since the observed supersonic line widths in giant molecular clouds (GMCs) are due 
to disordered macroscopic motions that contribute to dynamical support rather than to large-scale 
infall ( Zuckerman Evans 1974| ). While the large scale turbulent motions within GMCs lead 



to evolution of the shape of the cloud ( Ballesteros-Paredes et al 1999| ), they do not lead to an 



overall expansion or contraction on a dynamical time scale. Correspondingly, the lifetime of a 
GMC typically exceeds its free-fall timescale (~ 10 6 ' 5 y) by about an order of magnitude ( piitz fc 



Shu 1980] ). To a good approximation, then, most observed GMCs are in a steady state. Structures 



within GMCs also appear to be in approximate equilibrium: The density in and around a typical 
ammonia core within a GMC varies roug hly as r~ 2 (e.g., |Snell 198U [Fulkerson fc Clark 1984 



Zhou et al. 1990| ; [Ladd et al. 1991 ), a relationship that is consistent with a simple, stably stratified 



model. Turner (1993) finds that the majority of cores he observes in high-latitude cirrus clouds 
have characteristics that are consistent with hydrostatic equilibrium. 

We find that spherical models accurately represent many of the properties of these equilibrium 
structures, aside, of course, from details about the multi-dimensional cloud geometry. In particular, 
the time-honored problem of the structure and gravitational stability of gas spheres bears upon 
several unresolved issues in star formation. What is the dependence of density and line width on size 
scale for a cloud in stable hydrostatic equilibrium? What are the initial conditions appropriate for 
gravitational collapse, particularly on large scales where thermal pressure is relatively unimportant? 
How can GMCs be stable at the high mean pressures often observed? 



1.1. Molecular Cloud Structure 

Because we must specify an equilibrium structure before performing a stability analysis, our 
first step is to determine the density and line width profiles of a model molecular cloud. For a cloud 
in hydrostatic equilibrium, these profiles can be determined theoretically once the nature of the 
pressure supporting the cloud is specified. We consider three sources of pressure: thermal motions, 
static magnetic fields, and turbulence. 

The thermal sound speed in a given molecular cloud is roughly constant and corresponds to 
a temperature of 10-30 K. Thermal motions typically have significant dynamical influence only 



in cloud cores ( Myers fc Benson 1983 ). This pressure component may dominate support in these 
limited regions if the static magnetic field is nearly straight and uniform or if the field strength is 
low. 

Typically, though, static magnetic fields play a substantial role in stabilizing molecular clouds 
(Mouschovias 1976b; Heiles et al. 1993). Magnetic support is particularly significant in regions that 



are forming clusters of stars. Since the mass of a star-forming core exceeds its thermal Jeans mass, 
a clump containing multiple cores must contain many thermal Jeans masses; the material in the 
clump thus owes much of its dynamical support to the static magnetic field. 

Turbulent motions provide another source of nonthermal support of molecular clouds. On large 
scales, turbulent motions have an energy density comparable to that of the static magnetic field 
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( Myers fc Goodman 1988|). These motions are supersonic, and are observed to be an increasing 
function of map size ([Larson 198"H |Cernicharo et al. 1985| ; |Myers fc Fuller 1992|). It was suggested 
some time ago that these motions are large-amplitude MHD waves ( Arons &: Max 1975| ); recent 
obser vations ([Myers &: Goodman 1988| ) and theoretical results ( ponazzola ct al. 1987| ; pudritz 



1990| ; |McKee &: Zweibel 1995 ) are consistent with this idea. Alfven waves, being non-compressive, 
are expected to be the dominant form of MHD waves present ( McKee fc Zweibel 1995 ). The 



pressure due to the waves includes both the dynamic pressure associated with the motion of the 
gas and the magnetic pressure associated with time-dependent magnetic field. The importance of 
the wave pressure increases with size scale for two reasons: first, there is an increase in the range of 
wavenumbers for the MHD waves ( ponazzola et al. 1987 ), and second, the density decreases with 
size ( [Larson 1981 ), which leads to an increase in the amplitude of the waves flFatuzzo fc Adams 
1993| ; |McKee &: Zweibel 1995 ). Numerical simulations confirm that MHD waves are an important 
contributor to the support of molecular clouds, but they also indicate that the waves damp out 
very rapidly ( pammie fe Ostriker 199*6 , Vazquez-Semadeni et al 1999 ; but see McKee 1999] ). 



The primary pressure components that support molecular clouds against gravitational collapse 
have been identified, and there is substantial agreement as to when each plays an important role: 
Static fields are important throughout molecular clouds, whereas thermal pressure is important 
primarily in small, dense regions and turbulent, or wave, pressure is important on large scales. 
However, there is no quantitative understanding as to how these three sources of pressure interact 
to produce the density and velocity structure that is observed in molecular clouds. To date there 
has been only one study of the case in which all three components are present ( |Lizano &: Shu 1989[) . 
In this work, as well as in subsequent work (Gehman et al, 1996; McLaughlin & Pudritz 1996), a 
phenomenological model is adopted for the wave pressure, in which it varies as the logarithm of the 
density (a "logatrope" ) . The disadvantage of such phenomenological models is that it is impossible 
to study the stability of the clouds in a physically motivated way. Our intention is to develop a 
theoretical framework in which the distinct pressure components can be modeled in such a manner 
that both the structure and the stability of the clouds can be determined. 



1.2. Precollapse Conditions in Star-Forming Regions 

The state of a cloud at the onset of collapse determines the development of the density and 
velocity profiles once collapse begins, the fragmentation of the cloud, the luminosity of the protostars 
that form, and the final stellar configuration. We can illustrate the effects of the precollapse 
conditions by comparing self-similar solutions for the collapse of an isothermal sphere, which is 
often used to model low-mass star formation. There is an infinite number of such solutions (Hunter 
1977), with the solution discovered by Larson (1969) and Penston (1969) and the solution found 
by Shu (1977) demonstrating the ranges of velocity profiles and densities that are possible. The 
Larson-Penston (LP) solution starts from a static, uniform gas that fills a large spherical volume. 
As the gas collapses under the influence of gravity, a rarefaction front propagates inward from the 
outer boundary at the sound speed. When the rarefaction reaches the origin, a condensed object 
forms. At this instant, the inflow is spatially constant at 3.3 times the sound speed; the density 
has the same scaling as a singular isothermal sphere (SIS), p oc r -2 , but its magnitude is 4.4 times 
greater. Shu (1977) considered the other extreme, in which collapse begins from the equilibrium 
state described by the SIS. This solution is "inside-out" in the sense that after a perturbation at 
the origin, information that internal regions are evolving dynamically propagates outward in an 
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expansion wave. Since the initial state is static, with p oc r -2 , core formation begins immediately, 
before a flow field develops. Following the formation of a condensed object in both the LP solution 
and the Shu solution, the region around this object acquires power-law profiles of the form p oc r -3 / 2 
and v oc r _1//2 ; however, the density and velocity at comparable points remain larger in the LP 
solution. As a result, the mass accretion rate and the accretion luminosity are larger in the LP 
solution than in the Shu solution. 

Results from studies of the dynamical evolution of collapsing isothermal clouds often resemble 
the analytic solutions described above. An isothermal sphere in an unstable, but non-singular, 
equilibrium has a density profile that is approximately flat (like the initial state of the LP solution) 
in a central region, but that resembles the r~ 2 form of the SIS in an outer region. As a result, 
simulations of collapse from such initial conditions ( [Hunter 1977 ; Fiedler & Mouschovias 1993 



[Foster fc Chevalier 1993| ) display densities and flows near their incipient condensed objects that 
closely resemble those in the LP solution. By contrast, the material in the outer part of the sphere 
maintains acoustic contact while adjusting to imbalances between pressure gradients and gravity; 
thus, in this region, the simulation more closely reproduces the zero flow beyond the expansion wave 
in the SIS collapse solution. If the part of the sphere that approximates p oc r -2 prior to collapse is 
sufficiently large, the accretion rate eventually approaches the constant mass accretion rate of the 
SIS (Foster & Chevalier 1993). Similar results are found for the gravitational collapse of a cloud 



supported by a magnetic field undergoing ambipolar diffusion (Saner, McKee, & Stahler 1997). 
The initial state of the cloud also affects the stability of the collapse with respect to fragmentation, 



an essential process in the formation of star systems and clusters ( Boss 1988 ). The collapse of a 
uniform-density sphere is unstable to fragmentation, while that of an isothermal sphere is much 
less so (Shu 1977; [Silk fc Suto 1988| ). Simulations {e.g., Boss 1987) confirm that a more centrally 



concentrated initial density profile inhibits fragmentation during dynamical collapse. 

This discussion illustrates the impact of the precollapse conditions on star formation in the 
idealized case of an isothermal sphere. In reality, the gas may not be isothermal, and magnetic 
fields and turbulent pressure both contribute to the support of the cloud. One simplification 
remains possible for low mass star formation, however: it is believed to proceed quiescently, as 



ambipolar diffusion gradually reduces the pressure support due to static magnetic fields ( Mestel 



Spitzer 1956[ ; [Shu et al. 19871 ; |Mouschovias 199 1| ). In this case, the collapse proceeds from a critical 



point — an equilibrium state on the verge of gravitational instability. By this criterion, neither a 
homogeneous sphere nor an SIS are accurate representations of the precollapse state, since the first 
is not an equilibrium and the second is an unstable equilibrium. Nevertheless, highly centrally 
concentrated initial states are possible; calculations indicate that static magnetic fields provide a 
substantial stabilizing influence both in terms of the critical mass they support and the maximum 
center-to-surface density ratio that they allow ( |Mouschovias 1976b| ; Tomisaka, Ikeuchi, & Nakamura 
1988b [hereafter TIN]; [Lizano & Shu 1989|) . One of the objectives of this paper is to develop a 



framework for determining how these equilibria, particularly those at a critical point, are affected 
by Alfven waves. 



1.3. Pressures in Molecular Clouds 



GMCs appear to maintain a mean pressure that is substantially above the pressure in the local 
interstellar medium. Observations indicate that the mean kinetic pressures of nearby GMCs range 
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from 4 x 10 4 to 2.4 x 10 5 K cm" 3 QBertoldi fc McKee 1992|) , and the typical GMC in the molecular 
ring of the Galaxy has a mean pressure of 3 x 10 5 K cm based on the column densities inferred by 
Solomon et al. (1987). By comparison, the pressure in the local interstellar medium is only about 
2 x 10 4 K cm~ 3 ( |Boulares fc Cox 199C ; we have subtracted out the cosmic ray pressure since it is 



approximately uniform and cannot contribute to the support of GMCs) 

The origin of the relatively high mean pressures in GMCs is not understood. Consider the 
simple, well-studied example of the isothermal sphere: The maximum ratio of mean-to-surface 
pressure for a stable isothermal sphere is only about 2.5 (e.g., Spitzcr 1968j ), well below what 



is deduced from observations. The polytropes suggested by Maloney (1988) have even smaller 
pressure ratios. In the magnetized models of TIN, the maximum stable pressure ratio is also about 
3. Center-to-surface pressure ratios can be larger, but are also often limited — for example, a stable 
isothermal sphere has a maximum pressure drop of 14. As we shall see in §|5L however, it is possible 
to achieve larger center-to-surface pressure drops (see also Curry & McKee 1999). 

The resolution of this issue is crucial. Because the mean pressure of GMCs is representative of 
the pressure of the star-forming regions within them, it is intimately coupled to the star-formation 
process. The mean pressure determines, in part, the mean density in the star-forming region, which 
in turn affects the density of the star cluster that forms from the cloud. The line width-size relation 
( Larson 1981| ) suggests that that the mean pressure of self-gravitating GMCs is constant: In GMCs, 



the line width-size relation has the form a oc R 1 / 2 , where a is the 1-D velocity dispersion and R 
the cloud radius (e.g., Solomon et al. 1987| ); for self-gravitating clouds, we also have a 2 oc M/R, so 



that M oc R 2 , and hence the mean pressure, P ~ pa 2 oc (M /R 3 )R oc R°, should be the same for all 
GMCs. Chieze (1987) and Elmegreen (1989) have suggested that this constancy can be understood 
if the ambient pressure of GMCs is constant since, in their models, P is proportional to the ambient 
pressure. 

For both problems — the origin of the line width-size relation and the stellar number density 
ultimately produced by the cloud — an understanding of the high mean pressure of GMCs is essential. 
Part of the explanation lies in the fact that the pressure in the molecular gas is increased by the 
weight of the overlying atomic gas. Elmegreen (1989) pointed out the importance of this effect for 
the atomic hydrogen; Holliman (1995) showed that the effect is even larger if one is attempting to 
explain CO data on molecular clouds, since the layer in which the carbon is atomic is considerably 
thicker (and therefore heavier) than the atomic hydrogen layer. However, it is not clear that this 
effect is sufficient. Two other effects that could contribute to the large pressures observed in GMCs 
are (1) the ambient pressures of GMCs could be several times greater than the local interstellar 
value, perhaps due to the effects of massive star formation; or (2) the structure of the clouds could 
permit a mean pressure in the molecular gas significantly above the values allowed in existing 
models. The formalism developed herein can be used to evaluate the effect of the overlying atomic 
gas and to determine the mean pressure in the molecular gas. 



1.4. Previous Models of Self-Gravitating Gas Clouds 

In this paper we appeal to the relative simplicity of polytropic spherical models, in which 
P(r) oc /) 7p (r), to determine the structure and stability of self-gravitating gas clouds. Some of the 
earliest work on this problem was carried out by Ebert (1955) and Bonnor (1956), who investigated 
the stability of pressure-bounded, gravitating spheres supported exclusively by isothermal gas 
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pressure. These models become susceptible to dynamic instability when self-gravity induces a 
center-to-surface density ratio exceeding 14. Non-isothermal polytropes (7 P 7^ 1) have been 
considered by Shu et al (1972), Viala and Horedt (1974), and Chieze (1987) under the assumption 
that there is no internal heat flow (a locally adiabatic system; see §2.2). These results show that 
the maximum stable center-to-surface pressure ratio increases with 7 p . Maloney (1988) considered 
polytropes with 7 P < 1 in order to account for the observed line width-size size relation, since such 
polytropes have line widths that increase outwards. He assumed that some unspecified heating 
mechanism was able to maintain a constant central temperature when the cloud was compressed; 
with this assumption, clouds are stable to gravitational collapse. Lizano & Shu (1989), Gehman et 
al (1996), and McLaughlin &: Pudritz (1996) considered polytropes in the limit 7 P — > (logatropes) 
to model turbulent pressure in clouds. Isothermal polytropes (7 P = 1) with non-isothermal specific 
heats were analyzed by Yabushita (1968). 

Lynden-Bell and Wood (1968, hereafter LBW), used polytropes with 7 p = 1 to study self- 
gravitating star clouds (i.e., globular clusters) in order to understand the gravothermal catastrophe 
discovered by Antonov (1962). In contrast to the polytropic models of gas clouds cited above, 
internal heat flow is allowed, although there is no heat transfer from the cloud to its environment 
(a globally adiabatic system; see §2.2). Indeed, it is efficient thermal conduction inside the star 
cluster that establishes the isothermal temperature profile. They modeled the cluster as a gas 
with no internal degrees of freedom. If such a gas is confined by an ambient pressure, the onset 
of dynamical instability occurs at a center-to-surface density ratio of 389 (more precisely, we find 
389.6); the fact that the gas heats up when it contracts permits the cloud to evolve to much greater 
densities before it becomes unstable. However, for a cluster of stars it is appropriate to consider 
a gas confined to a fixed volume rather than by a fixed pressure since environment of the cluster 
performs no work on the stellar fluid. In this case, they showed that the maximum center-to-surface 
density for a stable configuration is 709. 

The structure of rotating, self-gravitating clouds has been discussed by Stahler (1983). Most 
clouds in equilibrium are thought to derive significant support from static magnetic fields, and 
magnetic braking should diminish the dynamical influence of rotation (Mestel & Paris 1984; 



Mouschovias 1987 ). The direct impact of static magnetic fields on the structure and stability of gas 
clouds is explored in the axisymmetric models of Mouschovias (1976a, b), Tomisaka et al. (1988a), 
TIN, and Lizano & Shu (1989). 

In §3 we present an approach for determining the structure of a spherical cloud supported by 
multiple pressure components. The stability analysis in §4 then allows us to determine the range 
of stable equilibria that can represent cloud structure. In many cases, these models are close to the 
critical equilibrium that defines the limit of gravitational stability, since it has been argued that 
ammonia cores (|Myers fc Benson 1983 ; Foster &: Chevalier 1993), massive star-forming clumps 



QBertoldi fe McKee 1992| ), and GMCs as a whole ( McKee 1989|) all verge on gravitational collapse 



2. Formulation of the Models 
2.1. Multi pressure Polytropes 



Our intention in this work is to determine the the structure and stability of gas clouds in 
hydrostatic equilibrium that are supported by several pressure components. The basic assumption 
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underlying our work is that such a quasi-static model can be used to model turbulent molecular 
clouds. We assume that in equilibrium, each pressure component satifies a polytropic relation of 
the form Pj(r) = -fTp^ 7 *" (r), where K p i is independent of position. In general, one can generalize 
the classical analysis of polytropes to multiple components in several ways: 



1. Multi-layered, or composite, polytropes have spatially distinct pressure components 
(Chandrasekhar 1939, p. 170), as in the core-envelope stellar models of Schdnberg & 
Chandrasekhar (1942). Applications to gas clouds are considered by Curry and McKee 
(1999). 

2. Multi-fluid polytropes have several different components that interact only gravitationally. 
Examples include multi-mass models for star clusters (Taff, Van Horn, Hansen, and Ross 
1975), models for clusters of galaxies including dark matter, and models for molecular clouds 
including embedded stars. 

3. Multi-pressure polytropes consist of a single self- gravitating fluid with several pressure 
components, so that the total pressure P(r) satisfies the relation 



It is this last type of polytrope that is of interest here. As discussed above, the pressure components 
relevant to molecular clouds are thermal pressure, static magnetic pressure, and the pressure due 
to MHD waves. 

We shall restrict our attention to spherical polytropes. Treating molecular clouds as spherical is 
clearly a substantial idealization, but it enables us to explore some of the essential physics underlying 
their structure. There are several effects that could lead to deviations from spherical symmetry: 
(1) Rotation. As discussed in §1.4 above, rotation is expected to have a relatively small effect on 
the structure of molecular clouds due to the effects of magnetic braking. (2) Tidal gravitational 
fields. These also appear to have a minor effect on Galactic molecular clouds (Scoville & Sanders 
1987). (3) Massive star formation. Massive stars can significantly disrupt molecular clouds, leading 
to violations of both hydrostatic equilibrium and spherical symmetry. Our models apply only to 
clouds that are not being disrupted by this process; fortunately, in the Galaxy there are many 
more GMCs than large OB associations (Williams & McKee 1997), so this will be a reasonable 
approximation for many clouds. (4) Anisotropic pressure. The stress due to an ordered magnetic 
field is intrinsically anisotropic, and necessarily leads to a violation of our assumption of spherical 
symmetry. Nonetheless, many of the qualitative features of magnetic stresses can be captured by 
treating them as spherically symmetric (Safier et al 1997), and it is possible to obtain reasonably 



accurate quantitative results for the structure of magnetized clouds using the approximation of 
a spherically symmetric magneti c stress flHolliman 1995|). (5) L arge scale turbulent motions can 



make clouds quite non-spherical ( Ballesteros-Paredes et al 1999|) ; however, the time average of the 



cloud shape is much closer to being spherical. 

Polytropes with a single pressure component are often described in terms of a parameter n 
such that p(r) oc T n (r); the polytropic index 7 p is related to n by 

7 P = 1 + -• (2) 

n 
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"Negative-index polytropes" with n < have been discussed by Shu et al (1972), Viala & Horedt 
(1974), and Maloney (1988). Such polytropes are hotter on the outside than at the center, just the 
opposite of a star. An example of such a polytrope that is relevant for our models is the pressure 
due to Alfven waves, which satisfy 7 P = 1/2 (McKee & Zweibel 1995), corresponding to n = —2. 
We shall assume that 7 P > 0, since polytropes with 7 p < are unstable for all values of the external 
pressure (Viala & Horedt 1974). 



2.2. Locally Adiabatic and Globally Adiabatic Pressure Components 

In order to determine the stability of multi-pressure polytropes, we assume that the response 
of each pressure component to an adiabatic perturbation is the same as that of an ideal gas that is 
either "locally adiabatic" or "globally adiabatic." Physically, a locally adiabatic pressure component 
has a time scale for internal heat transfer that is long compared to the dynamical time scale, whereas 
a globally adiabatic component has a short time scale for internal heat transfer. The distinction 
between local and global adiabaticity is necessitated by the fact that we are treating self-gravitating 
systems that are inhomogeneous. 

The response of a locally adiabatic pressure component to an adiabatic density perturbation 
satisfies the usual adiabatic relation, 

6 In Pi=jiS In p. (3) 
The specific entropy associated with a pressure component in our ideal gas model is 

k , f Pi\ 

Si = In — - + const (4) 



Ji-l VP 7 * 

(e.g., Landau & Lifschitz 1958). The entropy is thus determined by the entropy parameter 

Ki = 4; (5) 

pit 

this does not change with time for a locally adiabatic pressure component. Just as in the analysis of 
the stability of stars flLedoux 1965| ), the adiabatic indexes ji that describe the temporal variations of 
the pressure components need not be the same as the polytropic indexes 7 P j that describe the spatial 
variations. The polytropic indexes 7 P j help to determine the equilibrium structure of a cloud (§3); 
stability depends additionally on the adiabatic indexes ji (§4). Pressure components with ji = 7 p j 
are isentropic, since they have a spatially constant specific entropy that remains constant during 
an adiabatic perturbation: Ki = Pi(r) p(r)~ li oc ^(r) 71 "" 7 * is spatially and temporally constant 
for 7j = 7 p j. Chandrasekhar (1939) termed such components as being in convective adiabatic 
equilibrium since efficient convection produces such a structure in a thermally insulated system. 
On the other hand, if a locally adiabatic component has ji ^ 7 p j, then it cannot remain polytropic 
after a pressure perturbation; Yabushita (1986) has considered this case for spatially isothermal 
spheres (7 P = 1). 

If internal heat transfer is significant for a pressure component, then we assume that the 
component is globally adiabatic. Such a pressure component retains its polytropic form during 
an adiabatic perturbation due to efficient internal heat flow in the cloud. The distinction between 
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locally adiabatic and globally adiabatic components is important only for 7$ ^ 7 p j; if the component 
is isentropic, then no heat transfer is required to maintain the polytropic structure, and the 
component is locally adiabatic as well. The model of a globular cluster considered by LBW, in 
which the stars interact dynamically so as to maintain an isothermal distribution, provides an 



example of a globally adiabatic system. In our work, Alfven waves, which have 7 P « = 1/2 ( Walen 



1944j |Weinberg 1962]) and 7, = 3/2 flMcKee fc Zweibel 1995|) , are modeled as a globally adiabatic 
pressure component. (Note that since 7 P j 7^ 1, the "heat flow" associated with the Alfven waves 
does not make the gas isothermal.) The polytropes considered by Maloney (1988), like our model 
of a polytrope supported by Alfven waves, have 7 P j < 1 and are not locally adiabatic; however, 
Maloney's models require an injection of heat to maintain a constant central temperature during a 
compression, whereas our model redistributes heat that is already present. 

It is important to distinguish the physical pressure components from the model pressure 
components. We model each pressure component as a thermally insulated, ideal gas with a local 
ratio of specific heats equal to 7$. If the pressure component is actually thermally insulated, then 
7i is the ratio of specific heats, as usual. On the other hand, if the gas is subject to heating and 
cooling, then it can be modeled as a locally adiabatic component with an adiabatic index 7, that 
is generally not equal to the physical ratio of specific heats; the effects of heating and cooling are 
treated by invoking hypothetical internal degrees of freedom. For example, if the rate for heating 
is proportional to oc nT a and the rate for cooling varies as oc n 2 T b , then it can be shown that, in 
equilibrium, the model has 7^ = 1 + l/(a — b), whereas the gas (if it is monatomic) actually has a 
ratio of specific heats of 5/3. In such a case, Pi(r,t) depends only on density, and as a result, the 
pressure component is isentropic (7, = 7 P i). 



3. Structure Equations 

The elegance of polytropic models of self-gravitating clouds is due to the small number of 
physical constants and parameters that are needed to specify the structure of the cloud. The 
structure is determined by the mass equation, 

— = Anr 2 p, (6) 

and the equation of hydrostatic equilibrium, 

dP GM 



P- (7) 



dr r 2 

Note that these equations remain valid even if the sphere evolves quasistatically. 

It is customary to combine these equations into a single second order differential equation, the 
Lane-Emden equation (e.g., Chandrasekhar 1939| ; see Appendix [A|) . However, this equation admits 



a homology transformation, so that one of the two boundary conditions serves merely to set the 
density scale ( |Chandrasekhar 1939 ). As a result, it is possible to write this equation as a first order 



equation in terms of the scale-free variables 



4vrpr 3 . , GM(r)p 
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with a separate first order equation for the dimensionless radius ( Chandrasekhar 1939j ). 



We have chosen a variation of this approach, in which we place the two first order equations 
on an equal footing. We adopt 

X = HPc/p) (9) 

as an independent variable, where p c is the central density of the equilibrium cloud (i.e., it does 
not vary when the cloud is perturbed), x is generally a monotonically increasing function of r: 
Since P(r) is monotonically decreasing in a self-gravitating cloud, p{r) will also be monotonically 
decreasing provided only that the pressure varies as a positive power of the density (7 p j > 0), which 
is necessary for stability in any case (Viala and Horedt 1974). 

As dependent variables, we adopt one that is proportional to the mass measured in Jeans 
masses (cf Stahler 1983), 

M(r) 

^ = c 3 (r)/[G 3 p(r)]V2' (10) 
and one proportional to the radius measured in Jeans lengths, 

T 

A ^c(r)/[Gp(r)]V2' < U ) 

where c = {P/p) 1 / 2 is the isothermal sound speed. More specifically, p ~ M(r)/Mj(r), where 
Mj(r) is the local generalized Jeans mass — i.e., the maximum mass that the concerted action of all 
pressure components can maintain in equilibrium, given the values of specified parameters, such as 
the entropy and ambient pressure. Similarly, A ~ r/Rj(r), where Rj(r) is the radius of a uniform- 
density sphere containing Mj(r). These variables are related to the standard homology variables 
by 

u= , (n + l)v = ^. (12) 

p A 

An advantage of our dependent variables (and of the homology variables) is that they depend upon 
the properties at the surface of the cloud, not at the center. For interstellar clouds, the surface 
pressure is generally known, whereas the conditions at the center of the cloud are not. Further 
discussion of the relationship between our variables and the standard nondimensional variables can 
be found in Appendix |A|. 

Since individual polytropic pressure components can be expressed as Pi = K p i[p c ex.p(— x)] 7p S 
the total pressure is 

P = Y / K piP 2 pi e-^. (13) 

i 

Since din p = —dx, the overall polytropic index is 

_ d\nP _ dlnP_ i dine 2 
^ p d In p dx dx 

For a single pressure component, 7 p is spatially constant for polytropes. However, for a multi- 



component system, 7 p can be a function of position, since equations ( |i~3P and (14) imply that 

EK piP 2 pi e-^x ^\P 



= ^ b .^ P i - 7n ,.v = L ? > (Polytrope), (15) 
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where we have indicated that this equation is valid only for polytropes, which have K p i constant. 
In terms of our dimensionless variables, the mass equation (||) yields 

dlnM A 4 , , 

— — = 47r 7p — , 16 

whereas equation (|7j) becomes 

dlnr A 

—J— = Jp~- ( 17 ) 

Appearances of In M and In r in ( |i~6| ) and ( |i~7D can now be eliminated in favor of expressions involving 
p and A with the aid of equations ( |l0|) and (11), 

dkiM = din n + 3dln c - ^dlnp, (18) 

and 

dlnr = dlnX + cilnc — — din p. (19) 

We then use equations (|i~4|), (p^), and (|l~9l) in equations (|l^) and ( [l7|) to obtain the structure 
equations for multi-pressure polytropes: 

d In a A 4 1 , , . , . 

_Jf = fe7p ---(4-3 7p ), (20) 

and 

<i In A A 1 , , . 

-— = 7p --- 2- 7p . 21 
dx p 2 

If desired, one can obtain a single first order equation by taking the ratio of these two equations; 
once //(A), say, is found from integrating this equation, the density can be evaluated by integrating 
equation (^T|). 

The structure equations are integrated outward from the center of the cloud, which has 
/i = A = and % = 0. In order to interpret X/p near the center, we also need a relation between 
these two variables there. Since the cloud is not singular, the density approaches a constant at the 
center, and as a result 

u Mir) 47r 

4 = -V T (A 0) - (22) 

Equations (20) and (21), together with equation ( ]T5| ) for 7 p , describe the the general structure 
of multi-pressure polytropes. By inspecting these equations, we see that for a given 7 p , equilibria 
are distinguished only by how far integration is carried in x and, thus, by how centrally concentrated 
they are. Figure |] shows the structure of a spatially isothermal sphere; the origin corresponds to 
the center of the sphere, and the surface can be at any point on the curve. The stability of the 
sphere is discussed in §4 below. 
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3.1. Singular Polytropic Spheres 

Just as the singular isothermal sphere is very useful in studying the structure and dynamics 
of isothermal clouds (Shu 1977), so the singular polytropic sphere (which we label "SPS") is useful 
in studying the structure of polytropic clouds. The basic relations governing such spheres have 
been given by Chandrasekhar (1939) and, more recently, by Stahler (1983) and by McLaughlin 
and Pudritz (1996). Singular spheres have power-law profiles for the density, etc. The equation of 
hydrostatic equilibrium (^\) implies 

/90Cr -2/(2-7 P ) ) p oc r -2 7 p/(2-7 P ) ) c K r (l-7 P )/(2-7 P ) . ( 2 3) 

Observe that the sound speed c increases outward for 7p < 1, as discussed in §2. In order for the 
mass to be finite, we require p to fall off more slowly than 1/r 3 , which implies 7p < 4/3. One can 
then readily show that an SPS is described by single values of p and A 

/'< 
A 



(4-3 7p )^ 7p 



1/2-, 3/2 



2\V 2 x 

vrj (2- 7 p) 2 ' 

1 \V2 (4-3 7p )V2 7p i/2 



2tt 



7p 



(24) 
(25) 



An important parameter describing the thermal structure of a cloud is the ratio of the mean 
square value of c, 

(c 2 ) = T7 I <? dM = T7 I PdV = - ( 26 ) 
x ' MJ MJ p v ' 

to the surface value, where P and p are the volume-averaged pressure and density, respectively. 
For an SPS, we find 

_ (c 2 ) 4-3 7p 6 
^ = ^ = 6^' (7p< 5 } ' (27) 

where we have added the subscript "s" on c 2 to emphasize that it is measured at the surface of the 
cloud. 

Note that ip, which is equivalent to the mean energy per gram divided by the surface value, 
diverges as 7p — > 6/5: for 7p > 6/5, the energy of an SPS is concentrated at the center. Physical 
polytropes must satisfy the boundary condition dP/dr = at the origin, since the gravitational 
force vanishes there. The fact that SPSs do not satisfy this condition is not important for 7p < 6/5, 
since the mass and energy of such polytropes are dominated by the outer layers. As a result, the 
solutions for physical polytropes approach the SPS solutions for large values of x- Ori the other 
hand, for 7p > 6/5, the fact that SPSs do not satisfy the proper central boundary condition means 
that they cannot serve as approximations for physical polytropes. For example, equation (|2^) shows 
that an SPS with 7p = 6/5 has p oc r -5 / 2 , whereas the analytic solution for this case gives p oc r~ 5 
at large radii (Chandrasekhar 1939). 

The mean density in an SPS is 

p _ 3p 3(2 - 7p ) 4 
p s 4vrA' j (4 - 3 7p ) 3 
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whereas the mean pressure is 

P 3/# 3(2 - 7p ) 6 

(6 = fr£) (7p < 5 } - (29) 

(In each case, the first equation is general and the expression after — * is for the SPS.) As emphasized 
in the Introduction, polytropes do not have large mean density or pressure contrasts with the 
ambient medium: for 7 P < 1, both are less than a factor 3 for an SPS. 

The gravitational energy W can be described by the parameter a defined by 

3 GM 2 . . 

W = --a . (30) 

5 r 



For an SPS, we have 



:] " ;A %, d P <h- (3i) 



5 7p / 3 v/p 5' 



3.2. The Virial Theorem 



Before passing on to a consideration of the stability of multi-pressure polytropes, it is 
worthwhile to determine the implications of the virial theorem for our problem. The virial theorem 
for an unmagnetized gas sphere is 



3 / PdV - 3PV + W = 0. 



(32) 



Note that this equation applies at any point within the sphere as well as at its surface. In our 
notation, this becomes 

3 



3?/>A/u — 4tt\ air 

5 



0. 



(33) 



This can be solved to give an explicit relation for /u(A) in terms of two parameters, a and tp, that 
are usually of order unity: 

1ft \2\V2l 

lo7raA \ 



5ip\ 
~2a 



1± 1 



15^ 



(34) 



Note that [i has the correct limit ( p2| ) as A — > if the minus sign is chosen (since a and ip — > 1 as 
A — > 0). The solution with the minus sign is the correct solution for all values of A for 7 P < 0.82; 
it is correct for all stable isentropic spheres for 7 P < 0.90. For 7 P > 0.82, the solution can switch 
sign as A increases: For 0.87 > 7p > 0.82 there can be multiple sign changes, whereas for 7 p > 0.87 
there is only one change in sign, from minus to plus. Equation ( |34|) is exact, and it is not restricted 
to polytropes. 



4. The Stability of Multi pressure Polytropes 
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4.1. General analysis 

LBW have given a general discussion of the gravitational stability of spatially isothermal 
spheres (7 P = 1). They considered both the case of a sphere subject to a given pressure, which is 
appropriate for a cloud of gas, and that of a sphere confined to a fixed volume, which is a model for a 
cluster of stars. They focused on the cases j p = 7 = 1 (an isothermal sphere, which is the Bonnor- 
Ebert problem) and 7 P = 1, 7 = 5/3. Here, we shall generalize their discussion of pressure-bounded 
spheres to the case of arbitrary 7 and 7 P . Our objective is to find the criterion for determining 
the critical point, which is the most centrally concentrated structure that is gravitationally stable. 
Equivalently, the critical point gives the maximum mass for which an equilibrium exists at the 
specified conditions. 

The first law of thermodynamics for the cloud reads 5E = 5Q — P S SV, where E is the total 
energy of the cloud, 5Q is the net heat flow into the cloud, and P s is the pressure at the surface of 
the cloud. We shall focus on systems that are in equilibrium. Such systems are characterized by 
having a stationary value of the entropy subject to the given constraints. The equilibrium is stable 
if the entropy is a maximum. 

We wish to determine the stability criterion for a pressure-bounded gas cloud. The stability 
of a system is determined by its thermodynamic free energy. For an adiabatic, pressure-bounded 
system, this free energy is the enthalpy, H = E + P S V, the sum of the energy of the cloud and the 
work required to displace the volume occupied by the cloud. Note that P s is the pressure at the 
surface of the cloud, so that (in contrast to the energy) the enthalpy of the cloud is not calculated 
as an integral over the elements of gas within the cloud. Furthermore, as discussed in §2.2, an 
appropriate choice of the adiabatic index 7$ often allows one to model pressure components that 
are subject to heating and cooling as being adiabatic. Using the first law of thermodynamics, we 
find for the cloud 

5H = 5Q + V5P S . (35) 

Following LBW, we assess stability by considering a linear series of equilibria, an approach 
originally developed by Poincare (1885) for mechanical systems. Consider a series of gas spheres, all 
of which are in equilibrium and all of which represent clouds that have the same mass and entropy, 
but that are confined by an ambient pressure P s that varies monotonically along the series. (For 
example, this series of equilibria could result from the adiabatic compression or decompression of a 
cloud.) For spheres supported by more than one pressure component, we assume that the entropy 
of each component is constant along the series. These gas spheres need not be poly tropes; this 
analysis thus applies to perturbed locally adiabatic spheres, which are not polytropes if 7 7^ 7 P . In 
some cases, one will reach a critical point beyond which no equilibrium exists; i.e., 5P S = at the 
critical point. This point is an extremum and not a saddle point by construction, since we have 
assumed that there are no equilibria beyond the critical point; it will be a maximum or a minimum 
depending on whether the pressure must be increased or decreased to reach the critical point. Since 
the condition 5P S = cannot distinguish between an extremum and a saddle point, it is a necessary, 
but not sufficient, condition for a critical point. Now recall that in the series of equilibria that we 
have been considering no heat has entered or left the system {5Q = 0). It therefore follows from 
equation ([35]) that we also have 5H = at the critical point. 

Alternatively, one can fix the pressure and gradually change the entropy by allowing heat to 
flow into or out of the cloud. Again, depending on the values of 7 and 7 P , one may reach a critical 
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point beyond which no equilibrium exists. In this case 5Q = when one reaches the critical point; 
and since we assumed 5P S = 0, it follows that 5H = at the critical point for this sequence of 
equilibria as well. For simplicity, however, we shall focus on adiabatic pressure perturbations. Even 
though we assume that the perturbations in our model are adiabatic, we can treat heat transfer 
from the environment by allowing 7 to differ from the physical ratio of specific heats (see §2.2). 

A complementary approach is to consider a sequence of spheres with a monotonically varying 
mass, all with the same entropy (or entropy per unit mass) and all embedded in a medium of 
constant pressure. In this case, the critical point corresponds to the sphere with an extremum in 
the mass (5M = 0). For self-gravitating clouds, this extremum (if it exists) is the maximum stable 
mass under the specified conditions, and we label this M CT . Since 8M = at the critical point, it 
does not matter whether it is the entropy or entropy per unit mass that is held constant along the 
sequence.tl 

The existence of a critical point depends on the value of 7 and on the structure of the cloud, 
which is parameterized by 7 P for a polytrope. For example, isothermal spheres (7 = 7 P = 1; see 
Fig. ||), in which adiabatic compression increases the magnitude of the gravitational energy faster 
than the thermal energy, have a critical point, whereas spheres with 7 = 7 P > 4/3, in which the 
opposite is true, do not. Bonnor (1956) and Ebert (1955) used the condition 5P S = to determine 
the stability of isothermal spheres. 

Our approach in determining the stability of self-gravitating clouds is quite different from 
that of Maloney (1988) or McLaughlin and Pudritz (1996). Because these authors used a 
phenomenological equation of state, they could not define the entropy. As a result they were 
forced make an arbitrary assumption in specifying the applied perturbation; they chose to keep the 
central temperature constant. While this assumption may be plausible for the thermal pressure 
of cold molecular gas with T ~ 10 K, it is unjustified for wave pressure. Indeed, if the perturbed 
cloud is to remain polytropic for 7 P < 1, as they assumed, then 7 = 7 P and the cloud would cool 
as it is compressed; in order to maintain a constant central temperature, heat must be supplied 
throughout the cloud to counteract this adiabatic cooling. No mechanism has been advanced to 
account for such behavior, and we regard it as unphysical. 



4- 1.1. Isothermal Spheres 

The special case of an isothermal sphere is frequently considered in the interstellar literature. 
Such a sphere has a spatially constant temperature (7 P = 1) and its temperature remains constant 
if it is perturbed (7 = 1). By setting the ratio of specific heats for the gas to unity, we can 
treat isothermal pressure perturbations as being adiabatic; the work done on the gas goes into the 



1 The approach of consid ering a sequence of equilibria is also used in assessing the stability of degenerate stars 
( Shapiro & Teukolsky 1983). Stars have zero pressure at their surfaces, so the criterion 5P S = that we are using for 
clouds is irrelevant. Although degenerate stars have negligible entropy, they have a well-defined entropy parameter 
K and adiabatic index 7, which in general depend on the density but approach constants in the non-relativistic 
and extreme relativistic limits. The mass sequence can be arranged with either the central density or the radius as 
the independent parameter. The dependence of the equation of state on density leads to the possibility of multiple 
critical masses along the sequence (the Chandrasekhar mass and the maximum neutron star mass). In principle, such 
behavior can occur for multi-pressure polytropes as well. 
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hypothetical internal degrees of freedom. As a result, the critical point for this problem can be 
determined by extremizing H in the usual manner. It is customary, however, to use the Gibbs free 
energy G = H — TS and an adiabatic index equal to the actual ratio of specific heats. For an 
isothermal sphere, 5Q = T5S, so that equation ( |35| ) implies 5G = —S5T + V5P S . When 5T = 0, 
the critical point, which is characterized by 5P S = 0, also corresponds to the condition 5G = 0. 
Either way, one finds that the maximum possible density contrast for a stable isothermal sphere is 
14.0 (Bonnor 1956, Ebert 1955). 



4-1-2. Global Collapse vs. Core Collapse 

When an isothermal sphere collapses, the collapse is global: all flow is inward. On the other 
hand, when a globular cluster collapses, the collapse is restricted to the core; the envelope actually 
expands due to the transfer of heat from the core (LBW) . Such a heat transfer is generally required 
to effect core collapse, and as a result it cannot occur in clouds with only locally adiabatic pressure. 
For globally adiabatic clouds, the dividing line between core collapse and global collapse occurs 
at 7 = 4/3. Globally adiabatic spheres with 7 = 4/3 evolve homologously (Shapiro & Teukolsky 



1983 ): If one were to compress a 7 = 4/3 sphere, the rate at which the work done on the sphere 
is transformed into kinetic energy would exactly balance the rate of increase of the magnitude of 
the gravitational energy. Because the fraction of PdV work that is converted to kinetic energy 
increases with 7, compression of a sphere with 7 > 4/3 everywhere increases the ratio of kinetic to 
gravitational energy, reducing the center-to-surface density ratio p c /p s - Conversely, if the surface 
of such a globally adiabatic sphere were to expand, the ratio of kinetic to gravitational energy would 
decrease. If this evolution were carried sufficiently far, the cloud would become unstable, resulting 
in collapse of the core and dynamical expansion of the envelope. Thus, for 7 > 4/3, gravitational 
collapse is due to rarefaction of the cloud, not compression. 

Globally adiabatic clouds with 7 < 4/3 everywhere exhibit the opposite behavior. For such 
clouds, compression increases the magnitude of the gravitational energy faster than it does the 
kinetic energy, so that p c /p s increases. If the compression is sufficiently great, the cloud becomes 
unstable, and the entire cloud collapses. 



4.2. Equations for the Perturbed Structure 



In order to assess the stability of multi-pressure polytropes quantitatively, it is necessary to 
perturb the structure equations. This is done in Appendix [B|; we summarize the results here. For 
simplicity, we assume that there is only one globally adiabatic component. 



The equations governing the variations in p and A are 
dS In p 



and 



dSlnX 

dx 



dx 



<Hn7 p + 



2(51n7 p + (4-37p)(251nA-51n^), 



7 P 



+ 1 



5\np + 



37 P 



A 2 
p 3 



5 In A. 



(36) 



(37) 
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These equations involve the perturbation in the polytropic index, 5"f p , which is also evaluated 
in Appendix iBl. There we find it convenient to introduce a new adiabatic index 



E 



7i + 



(38) 



which is the weighted mean of the adiabatic indexes for the locally adiabatic components t and the 
polytropic index for the globally adiabatic component g. In terms of T, one finds 



Jin P 



3r 



win// ■[ ^ 



5]nK, 



py 



(39) 



This result together with equation ( B15 ) enable us to evaluate 7, the adiabatic index for a multi- 
pressure polytrope. In the simple case in which there is no globally adiabatic component, we 
find 

8lnP 6I11P K^fPi\ ,„ 

7^TT7T = -^— = r = E(^h. (^ = 0), (40) 



8lnp 



which is just the weighted mean of the individual values of 7$; thus, in this case, T reduces to 7. 
On the other hand, if there is a globally adiabatic pressure component but the locally adiabatic 
components are isentropic, then T = 7 P . 



In Appendix |B| we also define 



r' 



E 



liilp -7p») + 



P 



7 P9 (7p -7 PS )- 



(41) 



In terms of T and V , we find 



<51n7 p = - 



4( 7p -r)JlnA-2 7p-T + 



r 



4-3r 



<51n^ + 7 P g-7p 



3r' 



3r 



^ ) 



py ' 
(42) 

where 5 In K pg is the perturbation in the polytropic coefficient for the globally adiabatic component 
and is evaluated in Appendix O. If there are no globally adiabatic components, then this expression 
simplifies to 



5 In 7 P = - 

7 



4(7p-7)*lnA-2(7 P -7 + 



r 



4-37 



5 ln/i 



(43) 



The solution for the perturbed structure is then given by inserting equation (C17) into equation 
|), then inserting that into equations (|36| ) and ( |37[ ) for 8 In/j, and 8 In A, and then solving them 
simultaneously with equations ( p0|) and (21) for /i and A. This procedure is more complicated than 
the conventional treatment of the stability of hydrostatic clouds, which permits the stability to 
be determined from the solution of a single second order equation involving the spatial gradient 
of the adiabatic index 7 (Ledoux 1965). The additional complexity arises both because we have 
dropped the assumption that the perturbations are locally adiabatic, and because we have in effect 
evaluated the gradient of 7. 



As we found in §4.1, the critical point is given by the condition that <51nP s 
then implies that the critical point is determined by 



0. Equation <M 



2 

T 



p 



8 In K 



P9 ■ 



(44) 
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This condition can be used to determine the critical mass for clouds supported by locally adiabatic 
pressure components, for those supported by globally adiabatic pressure components, and for multi- 
pressure polytropes that may be supported by both types of pressure. 



5. Locally Adiabatic Polytropes 

In multi-pressure polytrope models of GMCs, the thermal gas pressure and the magnetic 
stresses are modeled as locally adiabatic pressure components: The entropy associated with the 
component remains constant in each mass element during a perturbation. Here we first determine 
how to characterize the average entropy associated with a locally adiabatic pressure component, and 
then use this to describe the critical mass. As we shall see, non-isentropic clouds are particularly 
interesting since they can have large central-to-surface pressure ratios. 



5.1. Entropy of Locally Adiabatic Pressure Components 

For locally adiabatic pressure components (labelled by £), the entropy parameter Kg = Pg/ p 11 
remains constant in each mass element during the evolution of the cloud (see eq. [| for the relation 
between Kg and sg, the specific entropy associated with pressure component £). In general, Kg 
depends on position unless the system is isentropic (7^ = j p g). 

Particularly for the case in which magnetic fields are important, it is convenient to introduce 
a spherically symmetric reference state of constant pressure Pg tTe { (Mouschovias 1976a). We shall 
generalize Mouschovias' concept of a reference state by allowing the density in the reference state 
p re f to depend on position in the cloud; we do not assume that the reference state is a polytrope. 
For a locally adiabatic component, the entropy parameter is the same for each mass element in the 
cloud as it is for that mass element in the reference state, 

Kg{M) = 4^ = Pl r\ • (45) 

In the isentropic case (7^ = j p g), the reference density is constant, as Mouschovias assumed: since 
the entropy parameter Kg is identical to the factor K p g that governs the polytropic structure, it is 
constant; with both Kg and Pg tTe { constant, it follows that p Te f is as well. 

For non-isentropic systems, we define the average entropy parameter by 

(Kz) 1 ^ = ±J K\ hl dM. (46) 



Introducing the power 1/7^ into the definition enables us to evaluate the integral for the reference 
state, 

M7i pi, 



\ K l) ~ — TIZ. — = -^ITi V 47 ) 



where dM = p Te fdV m { and V ve { is the volume in the reference state. Thus k\^ 1 is a measure of the 
local entropy per unit mass, and {Kg) 1 /* 11 is a measure of the mean entropy per unit mass of an 
entire cloud. 
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For an isothermal gas, we have P t h = pc 2 h with cth =const, corresponding to 7^ = 1. As a 
result, the entropy parameter is simply 

K th = (K th ) = 4. (48) 

Next, consider a cloud supported by a poloidal magnetic field. We treat the magnetic field as 
a gas with an adiabatic index jb = 4/3. This relation is exact for a tangled field. For a poloidal 
field, it is approximate, but it ensures that the magnetic flux <E> = J litrBdr is conserved during 
a homologous, spherical compression, since then B oc p 2 / 3 and rBdr oc r 2 p 2 / 3 (dr /r) =const. The 
parameter that measures the mean specific entropy for the field is 

p t/ 4 / 3 r2 p4 
(K R ) = ' rGf oc rcf ref (49) 

^ B; M 4/3 K M 4/3 ^ 

in terms of quantities measured in the reference state. Since the flux $ is the same in the cloud 
and in the reference state, we can rewrite this as 

{KB) ^ K ^73' (50) 

where B = Q/irR 2 is the mean magnetic field and p is the mean density. It follows that the 
parameter that measures the total entropy of the cloud depends only on the magnetic flux, 

M(K B f/ 4 oc $ 3/2 . (51) 

Adiabatic evolution of the cloud thus corresponds to evolution with constant magnetic flux. 



5.2. The Critical Mass 

As discussed in §4.1, the critical mass is the mass of a cloud at the critical point: For given 
values of the ambient pressure and component entropies, there is no nearby equilibrium state with 
a mass greater than M CT . In order to determine the value of the critical mass, we use the approach 
adopted in §4, based on a sequence of equilibria with constant mass and entropy but variable 
pressure. (It should be kept in mind that this sequence of equilibria is not a sequence of polytropes 
for non-isentropic clouds.) For simplicity, we consider a cloud supported by a single pressure 
component. We can express the mass of the cloud (eq. |H]) in terms of the entropy parameter at 
the cloud surface, Ki s , as 

M = -J^P S ^-^) K ^, (52) 

where we have introduced the subscript "s" to emphasize that the quantities are evaluated at the 
surface of the cloud. The critical point occurs at an extremum of the surface pressure, 5P S = 0. If 
we adiabatically perturb a cloud of fixed mass, equation (^) shows that the critical point occurs 
at 

6\n/j, = 0. (53) 



This is just the condition implied by equation (44) in the absence of globally adiabatic components. 
The critical mass is then 

3/2 2 
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where number that depends on je and j p £. 

For a non- isothermal gas, it is convenient to have an expression for the critical mass in terms 
of (c 2 ) = M" 1 f c 2 dM, which is often more readily observable. To accomplish this, we re-express 
equation ( |10| ) as 

fi'(c 2 ) 3 / 2 



M ^^W' <55> 

where 

A* 



recall that ^ = (c 2 )/c 2 (see eq. 27). At the critical point, we have fi' CT = fi CI /ipcr 2 ■ For an isothermal 
cloud, tp = 1 and p,' cr = /U cr . 

It is important to keep in mind that M cr is the maximum stable mass for given values of (Kg) 
(or, equivalently, Kgg) and P s ; it is the maximum stable mass for given values of of c s and P s (or, 
equivalently, p s ) only for isentropic spheres. In §5.3.1 below we show that p, is indeed a maximum at 
the critical point for isentropic spheres. For non-isentropic spheres, which are discussed in greater 



detail below (§5.4), p CT is less than the value for isentropic spheres. This can be seen graphically 
in Figure [l], which portrays all the equilibria for isothermal spheres: As 7^ rises above unity, our 
numerical results show that the critical point moves down and to the left, following this curve as 
the maximum stable value of x increases. (This result for isothermal spheres was first found by 
Yabushita 1968| .) How is it then possible to satisfy the critical point condition 5 In p = at a point 



other than the maximum in Figure [l]? The answer is that the perturbed sphere is not a polytrope 
for 7£ 7^ 7 and is therefore not represented in this figure. 

The expressions for the critical mass above depend on the specific entropy. It is also possible to 
express the critical mass in terms of the total entropy in the cloud, which is related to M(K() X I^ 1 . 
We define this alternative form for the critical mass, which we label M cr> t, by 

/M crt \ 3 Mcr , s 



V M ) M 

One can show that M cr t indeed depends on M and (Kg) only through the combination M (K^) 1 ^ 11 
with the aid of equation (|52|); it also depends on the numerical parameter (Kf s / {Ki)) cr . At the 
critical mass (M = M cr ), the two forms for the critical mass are identical (M cr = M cr t ). 

For an isothermal gas, the critical mass is the Bonnor-Ebert mass, 

M CT = M BE = 1-1822 f C ^ h 1/2 ) , (58) 

which follows from equation ( [54] ) with fi CT = 1.1822. The alternate form for the critical mass, M cr t , 

,2 

-tv 



depends on the total thermal energy in the cloud, Mc? h . 



Magnetically supported clouds have 7 = 4/3 as discussed above. The critical mass can thus 
be expressed as 

= Mcr a [WE) J cr ~CPI 2 ~ " G^V' (59) 
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where we have dropped the numerical factor (. . .) cr in the last step. Note that because 7 = 4/3, the 
critical mass (eq. |52] evaluated at the critical point) is independent of the pressure: for M > Mb, 
clouds are magnetically supercritical and cannot be prevented from collapsing by the magnetic 
field, whereas the contrary is true for M < Mb, under the assumption of flux freezing (cf. Shu et 
al 1987). The alternative form for the critical mass in terms of the parameter that measures the 
total entropy is 

M, = M cr>t = (M*M a )W oc ( M2 g;/ /2 ) V3 oc (60) 

The numerical value for M$ depends on the distribution of the flux in the cloud; for the standard 
case in which a uniform field threads a spherical cloud, it is given as 0.126$/G 1/2 by Mouschovias 
& Spitzer (1976). The relation between the two forms of the magnetic critical mass appears to 
have been first given by Mouschovias and Spitzer (1976); equation ( ]57|) shows that this relation is 
a general one that applies to all forms of support by an adiabatic pressure. 



5.3. Results for Single Component, Isentropic Polytropes (7 = 7 P ) 

The pressure components in interstellar clouds are "soft", with polytropic indexes 7 p j < 4/3, 
and this constrains the mass, density and pressure of self-gravitating interstellar clouds; for 
example, a stable isothermal sphere must be less than the Jeans mass and has a maximum density 
ratio of 14.04. Here we consider the limitations on the mass of single component, isentropic 
polytropic polytropes with 7 p < 4/3. 

The stability analysis of isentropic polytropes, previously considered by Shu et al (1987), Viala 
& Horedt (1974) and Chieze (1987), is a subset of our work, and we have plotted the results in 
Figures || and [|. Note that the critical density and pressure ratios increase with 7. The points 
marked on these figures denote the values for an isothermal sphere, the case originally considered 
by Bonnor (1956) and Ebert (1955). When 7 > 4/3, an isentropic polytrope is unconditionally 
stable. Such a polytrope can be stable even if the density and pressure vanish at its surface, and 
as such they can provide a simple model for stars. Systems with 7 = 4/3 evolve homologously 
in response to either a pressure perturbation or an entropy perturbation ( (Shapiro fc Teukolsky 



1983 ). The critical configuration for such a polytrope has a vanishing density at the surface and 
is neutrally stable against collapse. For example, white dwarf stars supported by ultrarelativistic 
degeneracy pressure (7 = 4/3) are neutrally stable, their stability being dependent on mass but 
not radius. Polytropes with 7 = 4/3 and with a finite density at the surface are unconditionally 
stable against collapse. 



The critical point condition 5 In /x = (eq. |53j) enables us to derive a simple relation between 
A and fj, at the critical point for single component isentropic polytropes. Such clouds have constant 
7 P , so integration of equations (20) and ( |2l| ) yields a unique solution. The perturbed cloud is also 



isentropic, so the perturbed cloud must lie on this solution as well. As a result, isentropic clouds 
have identical spatial and Lagrangian variations, so that <51n/i = <iln/i = at the critical point. 



(Recall from the discussion in §4.1 that this extremum is a maximum.) Equation (EG) then implies 



(61) 
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Chieze (1987) previously obtained this relation (his equation A18) when he considered the stability 
of single component polytropes, under the implicit assumption that 7 = 7p . The equivalence of the 
equations follows from the variable transformations in Appendix |A]. With the aid of the expression 
for /i obtained from the virial theorem (eq. [34]) it is possible to evaluate /x cr explicitly, 

P« = n ,„_xi/9 7 — 1 ; — (7 P = 7J- (62) 



The parameters a cr and V'cr are of order unity for 7 P < 1.3; as 7p — > 4/3, we find a cr — ► 2.5 and 
// cr — > 4.555, whereas ?/> cr oc (4 — 37 P ) _1//4 grows without bound. One can readily show that equation 
(|62| ) agrees with equation (^4|) for an SPS if the values of a and ip for an SPS are inserted into 
this expression. Note that fj, CT — > as 7 P — > 0: the critical mass shrinks with 7 P since it becomes 
impossible to maintain a substantial pressure gradient as 7 P approaches zero. 



5.4. Results for non isentropic, locally adiabatic polytropes (7 7^ 7p ) 

5. 4-1- Stellar vs. pressure-confined polytropes 

The behavior of locally adiabatic polytropes is summarized in Figure |5[ Since such polytropes 
adhere to the Schwarzschild criterion for stability against convection, those with 7 < 7p are 
convectively unstable; we do not consider them further. Locally adiabatic polytropes can be divided 
into two classes: 

Stellar polytropes (7 > 4/3, 7 P > 6/5) can serve as models for stable stars. Polytropes with 
7 > 4/3 are unconditionally stable (Ledoux 1965), since a locally adiabatic gas with 7 > 4/3 
becomes hotter at every point under compression and therefore cannot collapse under its own 
gravity. (Polytropes with 7 = 4/3 and zero pressure at their boundaries are neutrally stable.) 
Furthermore, pol ytropes with 7 P > 6 /5 can have a zero-pressure boundary at finite mass and 
radius, like stars ( |Chandrasekhar 1939[ ). (Recall that polytropes with 7 P > 6/5 can have all their 
energy concentrated in the core — § |3.1| .) As the surface pressure approaches zero ( x ~~ > 00 ) the 
dimensionless mass /i — * for 6/5 < 7 P < 4/3, 4.555 for 7 P = 4/3, and 00 for 7p > 4/3, as can 
be inferred from Chandrasekhar (1939). For fixed 7, there is no limit to the dimensional mass of 
a stellar polytrope; stars have limited masses because 7 depends on the mass, approaching 4/3 for 
massive, radiation-dominated stars or for degenerate stars that approach the Chandrasekhar limit. 

Pressure-confined polytropes (7 < 4/3 and/ or < 7 P < 6/5): First consider polytropes 
with 7 P < 6/5. All such polytropes must be confined by the pressure of an ambient medium 
(Chandrasekhar 1939). In the limit of a large center-to-surface density ratio {x ^ 1)) the pressure 
at the surface of such a polytrope approaches that of an SPS, 

_ 7p (4-3 7p )c 4 
s ~ 2^G(2-7 P 2 )V [M) 

(Polytropes with 7 p = 6/5 do not satisfy this relation because their mass and energy are 
concentrated at the center, but they too must be pressure-confined if they have a finite radius.) 
Next, consider polytropes with 4/3 > 7p > 6/5 and 7 < 4/3. Such polytropes can have arbitrarily 
large values of x, but they become unstable for x > Xcr', as a result, stable polytropes in this region 
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of parameter space must also be pressure-confined. The mass and energy of these polytropes 
is concentrated near the center when x becomes large. Finally, polytropes with 7 = 4/3 must be 
pressure-confined if they are to be stable (and not just neutrally stable) against small perturbations. 
Polytropic models for thermal pressure in interstellar gas and for magnetic pressure are in the 
pressure-confined regime, and as we shall see below, polytropic models for Alfven waves are 
equivalent to locally adiabatic polytropes in this regime. Insofar as the pressure in interstellar 
clouds can be represented as the sum of these pressure components, it follows that polytropes used 
to model interstellar clouds must be pressure-confined. 



5.4-2. Limitations on the mass and radius of pressure-confined polytropes 

The parameters fj, CI and fj,' CI , which describe the critical mass of a cloud, are plotted in Figures 
H and [7|. Recall that for the isentropic case, fi CT is a local maximum in /i(x) f° r 7p < 4/3 (§ |5.3D ; 
numerical solutions show that this is a global maximum as well, so that /x cr (non — isentropic) < 
/i cr (isentropic). Furthermore, as indicated in Figure ||, /j, cr (isentropic) is a monotonically increasing 
function of 7 P , reaching /x cr = 4.555 at 7 P = 4/3. The maximum value of n' cr generally occurs for 
non-isentropic clouds, and so cannot be read off Figure 0. We find that it too is a monotonically 
increasing function of 7 P , reaching max(/x^ r ) = 1.686 at 7 P = 4/3. We conclude that all polytropes 
with 7 P < 4/3, whether isentropic or not, satisfy /x < 4.555 and // < 1.686, so that 



M < 4.555 



(G^s) 1 / 2 



M < 1.686 



*}3/2 



(G 3 Ps )V2 



(64) 



In terms of the surface pressure instead of surface density, we find 



M < 4.555 



(G3P S )V2 



, M < 1.409 



,2\2 



(<*) 



(G3P S )V2 



(65) 



If attention is restricted to negative-index polytropes, the same line of reasoning that led to an 
upper limit of 4.555 on \i for 7 P < 4/3 leads to an upper limit /i < /i cr < /«cr(7 = 7p = 1) = 1.1822 
for 7 P < 1. Our results show that for such polytropes // < m&x([i' CT ) < 1.1822 as well. Note that 
since /i cr decreases below ^ cr (isentropic) as 7 increases above 7 P , it is possible for /x to exceed /i cr 
for stable clouds (similarly, // can exceed fj! CT \ but this is by less than a factor 1.5 for 7 P < 1. 

Our results also set limits on the size of polytropes. As indicated in Figure [l], the maximum 
value of A occurs prior to the critical point. This maximum value increases monotonically with 7 p 
for 7p < 4/3, reaching 0.6603 at 7 P = 4/3. Therefore, all polytropes with 7 P < 4/3 satisfy 



r < 0.6603 



(Gps) 1 / 2 



0.6603 



(GP B )V2 



(66) 



For negative index polytropes, the upper limit on the radius has the same form, but the coefficient 
is reduced to 0.5142. 
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5.4-3. Pressure and density drops in non-isentropic polytropes 

As discussed in the Introduction, an important question to be addressed by the polytropic 
models is whether they can account for the observed densities and pressures in molecular clouds. 
First consider the mean densities and pressures. These quantities have been plotted for isentropic 
clouds in Figure |I[ and they diverge as 7 P — > 4/3. For 6/5 < 7 P < 4/3, the mean-to-surface 
density and pressure ratios can become arbitrarily large as 7 — > 4/3, since such polytropes have 
their mass and energy concentrated near the center, as discussed above (§ |5,4.1| ). For 7 P < 6/5, the 
mean-to-surface density and pressure ratios can no longer diverge. For negative index polytropes 
(7 P < 1), it is convenient to express our results for the mean densities and pressures in terms of 



the values for SPSs (eqs. pq, 29), 



i <- ^Mr < 67> 

We conclude that negative-index polytropes cannot have large mean-to-surface density or pressure 
ratios, regardless of the value of the adiabatic index 7. 

Next, consider the central-to surface density and pressure ratios, which can become quite large 
for non-isentropic clouds. For fixed 7 P , the maximum value of x f° r a stable equilibrium (i.e., Xcr) 
increases as 7 increases above 7 P , and it is possible for Xcr to reach infinity if 7 exceeds a value 
we label 700 . We can determine the value of 700 analytically by noting that the cloud becomes a 
singular polytropic sphere with constant values of fi and A as 7 — > 700 . After a little algebra, we 
find that the requirement that the critical point (5 In fi — > 0) occur in the limit x — > 00 implies 

_ 32 7p (2- 7 p) 6 
7°°- (6- 7 p) 2 (7p -5 j - m 

For 7 p = 1, this yields 700 = 32/25 = 1.28. Yabushita (1968) previously pointed out that spatially 
isothermal spheres with \ = 00 are unstable for 7 < 32/25. For 7 p = 6/5, 700 = 4/3; for 
4/3 > 7 P > 6/5, 700 remains at 4/3. Note that polytropes with 7 > 700 and 7 P < 6/5 satisfy 
equation ( p3| ) for the surface pressure, so the infinite pressure drop is attained only at infinite 
radius. Since actual clouds are finite in extent, they are pressure-confined as discussed above. 

There is a subtlety in the behavior of non-isentropic clouds with large density drops, however. 
Recall that as x ~^ Xcr, the pressure perturbation <51nP s — > according to the condition for the 
critical point (§ |4.1| ). For 7 p = 1.2, for example, the pressure perturbation within the sphere remains 
comparable to the central value 5 In P c throughout most of the volume and approaches only near 
the edge. For negative-index polytropes, though, the pressure perturbation can drop to very small 
values well inside the sphere. Consider the example of a sphere with 7 p = 0.5 and 7 = 0.79 (slightly 
smaller than 700 = 0.793). We find that 5 In P drops to 0.15 In P c at \ ~ 4, far smaller than Xcr — 32. 
Thus, for any such polytrope with x > 4, the surface perturbations must be very small in order 
to prevent the perturbations at the center from becoming nonlinear. Although our linear analysis 
cannot determine what would occur if the central pressure perturbation became nonlinear, it is 
quite possible that the cloud would become unstable. If so, the stability of the cloud in this case 
could be likened to that of a pencil balanced on a flattened point — as x increases, the size of the 
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flattened region shrinks so that it becomes susceptible to smaller and smaller perturbations, and 
as a result, such a cloud would be unlikely to survive in a medium with large pressure fluctuations. 
As 7 increases, the value of % a t which the pressure fluctuation becomes substantially smaller than 
the central value grows. We find that a sufficient condition for density drops of at least 10 3 to 
occur in clouds with 8lnP s /5\nP c > 0.1 is 7 > l.^jp 1 ^ 2 . For j p < 0.8, clouds with 7 just above 
7oo do not satisfy this condition even though they are formally stable for arbitrarily large density 
drops. It is interesting to note that insofar as logatropes can be represented by polytropes with 
7 P — > 0, they do become stable for very large density drops provided 7 is sufficiently large; the 
central temperature is not constant unless 7 = 1, however. 

Finally, we consider the response of locally adiabatic polytropes to pressure perturbations. 
For isothermal polytropes, M cr is reduced somewhat by an increase in the external pressure, 
M cr oc P s *' 2 (eq. |54|). For negative-index polytropes, however, M cr can decrease much more 



sharply with P s due to the decrease in c s (|Shu et al 1972| , |Tohline et al 1987|) . Since the decrease 



in the temperature is bounded (it is difficult to cool below 10 K in a typical molecular cloud, 
for example), it is more convenient to express the critical mass in terms of quantities after the 
compression, 

where c s j is the final value of c s , etc. This result shows that the reduction in the critical mass due 
to cooling, which reduces c s , can be large, but the reduction due to the compression is limited by 

— 1/2 

the weak P s * dependence. For example, in a radiative shock the final pressure is related to the 

initial value P s ,i by P s j = (fshock/c s ,i) 2 -fs,j- In spherical implosions even higher compressions, and 
correspondingly greater reductions in M cr , are possible ( [Tohline et al 1987 ), although in practice it 



may be difficult to maintain the high degree of spherical symmetry required to achieve very large 
compressions. 



6. Globally Adiabatic Polytropes 

When the flow of energy is significant on a dynamical time scale, the assumption that the 
perturbation is locally adiabatic breaks down. The perturbation can be approximated as being 
globally adiabatic provided that no heat is supplied to or removed from the cloud during the 
perturbation. This is the approximation that we adopt for the Alfven waves we are using as a 
model for the turbulence in molecular clouds. Precisely because heat can flow during a perturbation, 
determining the condition that the perturbation be adiabatic is non-trivial, and it is this problem 
we now address. 



6.1. Entropy of a Globally Adiabatic Pressure Component 

For a globally adiabatic pressure component g, heat is assumed to flow within the system so 
as to maintain the polytropic condition P g = K pg p^ P3 , but there is no exchange of heat with the 
surroundings, nor with the other pressure components. As a result of the heat flow, K pg may 
change with time, although it is independent of position. In the actual globally adiabatic system 
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(as opposed to our model of the system), this heat flow is reversible. For example, if a non- uniform 
cloud with Alfven vwaves is adiabatically compressed and then decompressed back to its original 
size, the Alfven waves will be unchanged by the process, since the wave action is conserved (Dewar 
1970). This leads to a problem: in our model for the Alfven waves, heat flows across a temperature 
gradient — an irreversible process. As a result, the entropy as conventionally defined may not be 
conserved for globally adiabatic processes. 

To define an entropy for a globally adiabatic pressure component, we again introduce a 
spherically symmetric reference state that is related to the cloud by a reversible adiabatic transition 
(c.f. Mouschovias 1976a). In this reference state, the globally adiabatic pressure component has a 
uniform pressure Pg )T ef] since globally adiabatic components are assumed to always be polytropic, 
it follows that the density p g re f is uniform as well. The entropy parameter for this component in 
the reference state is 

K 9 M - (71) 
Pref 

If the cloud then returns to its initial state by a reversible adiabatic process, the entropy will be 
unchanged. Thus, the entropy of a globally adiabatic pressure component is characterized by K g ^. 
In our case, we are interested in self-gravitating clouds, so the transition to the reference state can 
be visualized as being due to slowly decreasing the value of the gravitational constant G to zero. 

In Appendix |C| we evaluate K g ^ e t by extending an argument due to McKee & Zweibel (1995). 
There we demonstrate that K gTe f = (K g ), where (K g ) is the mean entropy parameter for a globally 
adiabatic pressure component, defined as 



(K g )\^-^J =— K} lp9 ~~' a) dM. (72) 

Observe that since 

K g = Cexp[( 7g - T)8 g /k], (73) 

where C is a constant and s g is the specific entropy of a gas element (eq. |j), it follows that the 
entropy parameter for a globally adiabatic pressure component is given by a well-defined integral 
of the specific entropy. However, in contrast to the case for a locally adiabatic pressure component 
(§5.1), the specific entropy s g may not be constant during a globally adiabatic change, due to 
internal heat flows; it is only the integral in equation ((7^) that remains constant. 

For a spatially isothermal component, for which there is no heat flow across a temperature 
gradient, the entropy parameter is related to the mean specific entropy (s g ) = M _1 J s g dM in a 



simple way. Let e = 7 P9 — 1 be small. Inserting K g (s g ) from equation ( |73| ) into equation (|72), we 
have 



C 



(Kg)™-™ = — I exp 



7g-M s 



Jvg- lg) k 

Expanding both sides of this equation for small e, we find 

7 S -1 



dM. (74) 



ln{K g ) = (s g ) + const. (75) 

Since (K g ) is constant for a spatially isothermal, globally adiabatic system, it follows that {s g ) is 
too (c.f. LBW). However, if the system is neither spatially isothermal nor isentropic (i.e., 7 P9 ^ 1 
and 7 ps ^ 7 5 ), then only {K g ) is constant. 
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The polytropic constant K pg evolves during a globally adiabatic change as 

K%- 1 = (K g )T»- 1 (<i)' r '^» (76) 

from equation (|C9|). This can also be expressed in terms of the density by noting that (? g = 
Kpgp^B' 1 , so that 

K pg = (K g )(p^~ 1 )^ . (77) 

This equation allows one to follow the evolution of the pressure of a globally adiabatic component, 
P g = K pg p lps , in terms the evolution of the density structure of the cloud. 



6.2. The Critical Mass 

To determine the critical point for a cloud supported by a globally adiabatic pressure 
component, we express the mass as 

p 3/2 p(3 7 p 9 -4)/2 7p9 K 2/ lpg 

*=Ar <?» ■ (78> 

Using this equation to evaluate the effect of a pressure variation, we find that the critical point 
(SP S = 0) occurs at 

5lnp + \-^-JSlnK pg = 0, (79) 

which is just what we found in equation (f44|). Determination of the location of the critical point 
requires evaluation of 5 In K pg at constant entropy, which is carried out in Appendix |C]. 



6.3. Results for Globally Adiabatic Clouds 

There are two key differences between the stability of globally adiabatic clouds and that of 
locally adiabatic clouds: First, there is no convective instability, since it is assumed that the heat 
flow is rapid compared to fluid motions; as a result, it is possible to have stable systems with 
7 9 < 7 P9 , although we have not identified one in practice. Second, as discussed in §4.1.2, clouds 
with 7 9 > 4/3, which are stable if they are locally adiabatic, can be unstable to core collapse. In 
core collapse, instability is triggered when a perturbation causes a sufficiently large heat flow from 
the center of the cloud to the envelope. This is the case that is relevant to clouds supported by 
Alfven waves, which have j w = 3/2 (see §p\^). Globally adiabatic polytropes with 7 P > 6/5 and 



7 9 > 4/3 are not subject to core collapse, however: As shown in §5/4, polytropes with j p > 6/5 
can have arbitrarily large density contrasts, and our numerical results show that they are stable 
for 7 9 > 4/3. 

Results for the center-to-surface pressure and density ratios, P c /P s and p c /p s , are plotted in 
Figures ^ and ||, respectively. Note that for the cases of greatest interest, in which 7 S > 7 P9 , these 
ratios are determined primarily by the value of j pg . Comparing these results with those for locally 
adiabatic polytropes discussed in §Bl we find that globally adiabatic polytropes have smaller critical 
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density and pressure contrasts (provided 7 > 7 P , the only case of relevance for locally adiabatic 
polytropes). 

It is important to note that the structure of a cloud supported by a globally adiabatic pressure 
component is determined by the value of 7 P and is therefore identical to that of a cloud supported 
by a locally adiabatic pressure component with the same value of -y p . For a locally adiabatic cloud, 
the density drop at the critical point (p c /Ps)ci rises smoothly from the isentropic value (Fig. ||) to 
00 as j£ rises from 7 p £ to 7oo. Figure |9] shows that (p c /p s ) cr increases above the isentropic value as 
7 increases above 7 P for a globally adiabatic component as well. Hence, for any globally adiabatic 
cloud with 7 g > 7 P , it is always possible to find a value of 7^, which we label 7^(equiv), such that 
the critical point of the locally adiabatic cloud is identical to that of the globally adiabatic cloud. 
The value of 7^(equiv) for the equivalent locally adiabatic cloud is in the range 700 > 7^(equiv)> 7 P . 
We have evaluated 7^(equiv) for several cases: For 7 = 3/2 and 7 P = 0.5 (the Alfven wave case-see 
below), 7^(equiv)= 0.60564; for 7 = 3/2 and 7 P = 0.1, 7^(equiv)= 0.11985; and for 7 = 5/3 and 
7 P = 1.0 (the LBW case), 7^ (equiv)= 1.22042. It must be kept in mind that although it is possible 
to locate the critical point in a globally adiabatic polytrope with a locally adiabatic equivalent, 
the nature of the instability in the two cases is quite different: for j g > 4/3, the globally adiabatic 
polytrope is subject to core collapse, which is impossible for a locally adiabatic polytrope. 

How does this result carry over to the case of a cloud with multiple pressure components? 
Just as in the single component case, the structure of the cloud is determined by the values of 
7 P for each of the components, with the values of and j g determining which of the structures 
corresponds to the critical point. The critical point for a cloud supported by one or more locally 
adiabatic pressure components plus a globally adiabatic pressure component is therefore the same 
as that of a cloud supported by the same locally adiabatic components plus a locally adiabatic 
component with a suitable 7^(equiv). Numerical calculations for two component clouds show that 
the value of 7^(equiv) is not constant, but instead depends on the fraction of the pressure in the 
globally adiabatic component. This is reasonable, since although the structure equations do not 



depend on the values of ji and j g , the equations for the perturbed structure do (§4.2). In contrast 
to the single component case, the value of 7^(equiv) can exceed 700, although for Alfven waves it 
is generally £ 1. 



6.4. Clouds Supported by Alfven Waves 

We now focus on Alfven waves, which we are using to model the turbulence in molecular 
clouds. As shown by Dewar (1970), Alfven waves exert an isotropic pressure and can thus provide 
support parallel to the background magnetic field. For small amplitude waves, McKee & Zweibel 
(1995) showed that 7 W = 3/2 and 7 PjW = 1/2; this latter result was originally derived by Walen 
(1944). While Alfven waves are undoubtedly an important component of the MHD waves that 
contribute to the turbulent pressure in molecular clouds, there are several limitations to our model 
that should be borne in mind. First, we assume that the waves are adiabatic, whereas simulations 
suggest that the waves damp rapidly ( [Vazquez-Semadeni et al 1999| ). Second, following Dewar, we 
assume that the waves are in the WKB limit (short wavelength), whereas the observations suggest 
that the largest amplitudes are in waves with wavelengths comparable to the size of the cloud; as 
a result, we expect our models to be more accurate in reproducing the mean value of the wave 
pressure than the value of the wave pressure at the edge of the cloud, for example. Finally, whereas 
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simulations of ID MHD turbulence by Gammie & Ostriker (1996) confirmed that 7 PjW = 0.5 for 
small amplitudes, they indicate that 7 PjW becomes significantly smaller for large amplitudes. It is 
not clear that this result can be applied to GMCs, however: much of the reduction in 7 PjW they 
found was due to the fact that the gas became highly clumped at large amplitude, and the wave 
pressure was not much larger inside the clumps than outside. However, the question for GMCs 
is whether Alfven waves can maintain a pressure gradient in a clumpy medium, and this would 
require a much larger scale simulation. 

It is possible to infer 7 PjW from observation, but the existing results are not definitive. If 
observations show that the density scales as p oc r~ p and the ID nonthermal velocity dispersion 
scales as a n t oc r q , then P w oc pa 2 t implies 

7 P .w = 1 - — • (80) 
V 

Myers' (1985) reanalysis of Larson's (1981) heterogeneous data gave p = 1.2 and q = 0.3, 
corresponding to 7 PjW = 0.5. When Caselli & Myers (1995) focused on low mass cores, they found 
p = 1.1 and q = 0.53, corresponding to 7 PjW ~ 0; such cores are believed to be largely supported by 
static magnetic fields and thermal pressure, however. (It should also be noted that their value of p 
applies well outside the inner, isothermal core.) For massive cores, where the non-thermal motions 
are more significant, they found p = 1.6 and q = 0.21, corresponding to 7 P>W ~ 0.75. With the 
exception of low mass cores, then, it appears that 7 PjW ~ 0.5 is not inconsistent with the data. 

Adopting 7 PjW = 0.5 and 7 W = 1.5, we find that the maximum possible pressure ratio for 
a cloud supported by Alfven waves is only 4.15; the corresponding density ratio of 17.26 is 
somewhat greater than that for a critical isothermal cloud. The polytropic constant varies as 
K pg = {c 2 g ) 2 /{K g } = (K g ){p- 1 / 2 )- 2 (eqs. ^ and \n§ . The critical mass for a cloud supported by 



Alfven waves can be expressed in terms of the non-thermal velocity dispersion a n t since the Alfven 
wave pressure, including the pressure of the fluctuating fields, is 

3 

= pel = -pal, (81) 

(McKee & Zweibel 1995). For 7w = 3/2 and 7 PjW = 1/2, we find p, cx = 0.21334, Vcr = ((c 2 g )/c 2 ) cr = 
{{al t )/al t ) CT = 0.71342, and p' CT = 0.35405. The critical mass in this case is then 

a 2 la 2 ) 3 / 2 

M w = 0.3919 2*7 = 0.6504 x nt/ 1/0 . (82) 

GV 2 P l /2 G 3 / 2 pl /2 

This result for the critical mass is less than the estimate of McKee &: Zweibel (1992) based on the 
virial theorem; they found a coefficient of 2.17 in the second expression. This discrepancy is due to 
the fact that the maximum possible value of p in equilibrium (independent of stability) is governed 
by 7p> which does not enter the virial theorem argument. In inferring the maximum mass allowed 
by the virial theorem, one implicitly assumes that the three terms in equation ( |3"2| ) can be of similar 
magnitude; however, for values of 7 P significantly less than 1, the pressure becomes approximately 
constant in the cloud and this is not allowed. 



Although the results of this paper are primarily for single component polytropes, we shall 
present two results for multi-pressure polytropes. We find that the critical mass for a cloud 
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consisting of a cloud supported by isothermal gas pressure and Alfven waves is approximately 

M cr = (M^ E 3 + M2/ 3 ) 3 / 2 , (83) 

where the effective velocity dispersion cr c g is given by 

4f = C?h + (/^w,cr/^BE,cr) 2/3 (4)- (85) 

For Alfven waves with 7 p = 1/2, the effective velocity dispersion becomes 

^ff = c t 2 h + 0.67<<4); (86) 

equation (|||) is accurate to within 2% in this case. Thus Alfven waves are less effective at supporting 
clouds despite the fact that their pressure is 1.5 times larger than pa\ t due to the fluctuating 
magnetic fields. 



Second, in Figure 10 we show the density distribution in a multi-pressure polytrope that 
reproduces some of the features observed in molecular clouds. We have chosen the polytropic 
index of the magnetic pressure to be j p ,b = 1-0 in order to illustrate the dramatic effect of a 
non-isentropic pressure component (recall that 7# = 4/3). The magnetic pressure is comparable 
with the Alfven wave pressure, and the thermal pressure is 10% of the magnetic pressure. The 
case shown is a subcritical cloud, with a center-to-surface density contrast of 545 compared to the 
critical value of 1737; this density drop is far larger than that of the Bonnor-Ebert sphere shown 
for comparison, which has a density drop of 14. On the other hand, the mean-to-surface density 
contrast of 2.59 for the multi-pressure polytrope is only slightly greater than the value of 2.46 for 
the Bonnor-Ebert sphere. The inclusion of the wave pressure is destablizing: in the absence of 
wave pressure, a cloud supported by a magnetic pressure with 7 Pj b = 1.0 and a thermal pressure 
Pth = 0.1 Pb is stable for an arbitrarily large density drop. 



7. Summary 

The structure of molecular clouds has a major effect on the star formation that occurs within 
them. Polytropic models, in which the pressure varies as a power of the density (P oc /? 7p ), are 
useful in the study of the structure of such clouds just as they are for the study of stellar structure 
(Chandrasekhar 1939). Like stars, molecular clouds and some of the clumps within them are 
self-gravitating gas clouds, but there are two important differences: molecular clouds have a soft 
equation of state (for example, the gas is often approximately isothermal), so that they must be 
confined by an external pressure; and they are highly turbulent, which leads to a complex internal 
structure and a non-spherical shape, both of which are time dependent. In order to determine 
the characteristics of stable molecular clouds and the conditions under which molecular clouds 
become unstable to gravitational collapse, we have assumed that the pressure in these clouds 
can be represented by a sum of polytropic pressures — a multi-pressure polytrope. The individual 
terms represent thermal gas pressure, magnetic pressure and turbulent pressure. We model the 
turbulent pressure as being due to a superposition of Alfven waves, which have a polytropic index 
7 P = 1/2. This is an example of a "negative-index polytrope" (so called because the index n 
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in 7p = 1 + 1/n is negative), which has long been known to account for the outward increase in 
line widths, both thermal (Shu et al 1972) and non-thermal (Maloney 1988). Our models assume 
hydrostatic equilibrium and spherical symmetry, both of which become better approximations if a 
time average is taken. Rotation is neglected, which appears to be a reasonably good approximation 
(Goodman et al 1993). 

Several features of our analysis bear mention. First, since interstellar clouds are confined by the 
pressure of the ambient medium, we describe the cloud in terms of variables (fi = M[G 3 / 2 p 1//2 /c 3 ] 
and A = r[G 1 / 2 p 1 / 2 / c], where c = [P/p] 1 ^ 2 ) that depend only on local properties of the cloud, not 
on the properties at the center of the cloud. The variable fi is proportional to the ratio of the mass 
to the local value of the Jeans mass, whereas A is proportional to the ratio of the radius to the local 
value of the Jeans radius. The resulting description is simpler than the standard one (see Appendix 
[A] and Stahler 1983). Next, we assume that the perturbations are adiabatic. As is often done in the 
study of stellar structure, we allow the adiabatic index ji, which measures how the pressure of the 
ith component responds to a perturbation in density, to differ from the polytropic index 7 p j. The 
structure of the cloud depends on the values of 7 P j, whereas the stability depends on the values of 
7i as well. Third, for pressure components that are not isentropic (i.e., those for which 7, 7^ 7 P i), 
we consider two cases: For a locally adiabatic component, such as the thermal pressure or the 
magnetic pressure, the entropy of each mass element remains constant during a perturbation, but 
as a result the perturbed cloud is no longer a polytrope. On the other hand, for a globally adiabatic 
component, the entropy flows so as to maintain the polytropic relation between the pressure and 
density. This case is appropriate for Alfven waves under the assumptions that the perturbation 
is applied quasistatically and that there are no sources or sinks for the waves. We describe both 
locally and globally adiabatic pressure components in terms of the entropy parameter Ki = P/ p 1 * . 
For a thermal gas, K{ is just c 2 h ; for the magnetic field, the appropriate average of Ki is related to 
the magnetic flux; and for Alfven waves, it is related to the wave action. In Appendix O, we show 
how to define the mean entropy parameter for a globally adiabatic pressure component so that it 
remains constant during an adiabatic change. 

The essential feature of the three pressure components that support molecular clouds is that 
their equations of state are soft, with a limited ability to resist gravity: The thermal gas is 
approximately isothermal (7 P ~ 7 ~ 1) in regions of high extinction. The magnetic field must 
have 7 = 4/3 in order to ensure flux freezing in spherical symmetry (§ |5.1| ); stability against the 
interchange instability requires 7 P < 7 = 4/3. The turbulent pressure is modeled with Alfven 
waves, which contribute a significant fraction of the pressure in MHD turbulence, and these waves 
can be rigorously shown to obey 7 = 3/2 and 7 P = 1/2 in the limit of small amplitude, small 
wavelength, and negligible damping (McKee & Zweibel 1995). In contrast to the gas and the field, 
Alfven waves are globally adiabatic, and clouds supported by Alfven waves are subject to core 



collapse if the pressure drop inside the clouds is too large. As discussed in §0, Alfven waves are 
equivalent to a locally adiabatic pressure component with 7^(equiv) 1 insofar as determining the 
conditions for gravitational collapse. 

The results in this paper are primarily for polytropes with only one pressure component; 
results for multi-pressure polytropes are deferred to another paper. We extend previous work on 
polytropes with a locally adiabatic pressure to the non-isentropic case (7 7^ 7 P ). The parameter 
space for the stability of locally adiabatic polytropes is portrayed in Figure ||. Clouds with 7 < 7 P 
are unstable to convection and are not considered here. Locally adiabatic polytropes are classified 
either as stellar polytropes (7 > 4/3 and 7 P > 6/5) or as pressure-confined polytropes. Stellar 
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polytropes can serve as models for stars since they are not subject to gravitational collapse and 
they can have a zero pressure boundary at a finite radius. On the other hand, pressure-confined 
polytropes must have a finite pressure at their surfaces, either because 7 P < 6/5, or in order to be 
stable against gravitational collapse with 7 < 4/3. A polytrope supported by any of the pressure 
components in a molecular cloud must be pressure-confined: In contrast to the case of a star, the 
stability of a molecular cloud is determined by conditions at its surface. 

We introduced the concept of a globally adiabatic pressure component as an idealized model 
for treating the turbulent motions in molecular clouds. Molecular clouds are magnetized, and 
as a result the turbulent motions can be considered to be a superposition of MHD waves. The 
pressure associated with wave motions cannot be locally adiabatic since the waves can move from 
one part of a cloud to another under the influence of a perturbation. However, in the absence of 
sources and sinks, the total entropy of the waves remains constant — the system is globally adiabatic. 
A gravitationally bound cluster of stars is another example of a globally adiabatic system. We 
have determined the appropriate entropy parameter for a non-uniform, globally adiabatic pressure 
component; for the spatially isothermal case (7 P = 1), it is just the exponential of the usual entropy 
evaluated by LBW. We then generalized the standard treatment of the stability of polytropes to 
include a globally adiabatic component with an arbitrary value of 7 p . In contrast to locally adiabatic 
systems, globally adiabatic systems can be unstable to gravitational collapse even when 7 > 4/3; 
provided 7 P < 6/5, they undergo core collapse beyond the critical point (the "gravo-thermal 
catastrophe" of LBW). This is in contrast to the phenomenological model used by Maloney (1988) 
in which gravitational instability is impossible. Our results for the critical pressure and density 
ratios are shown in Figures || and pt these are smaller than for the corresponding locally adiabatic 
cases. Alfven waves, which we have used to model for the turbulent pressure, are particularly 



simple because their pressure is isotropic; as shown in equations ( |84T ) and (8£) and as discussed at 



the end of S.4, the non-thermal motions associated with such waves are somewhat less effective in 



supporting a cloud than are the thermal motions of an isothermal gas. 

The critical mass — i.e., the maximum mass that is stable against gravitational collapse for a 
given entropy distribution and ambient pressure — is evaluated in terms of the parameters // cr and 

m„ = = i4$2£. (87) 

The values of /z cr and /jl' ct for isentropic polytropes are plotted in Figures || and ^; the values of fi 
and // for singular polytropic spheres are also shown. The values of ii cr and n' cr depend primarily 
on the value of 7 P ; both are of order unity for 7 P ~ 1 and approach zero as 7 P — > 0. The maximum 
value of \i for any polytrope with 7 P < 4/3 is 4.555, the value of /i cr at 7 P = 4/3. Negative-index 
polytropes have /i < 1.1822, the value of // cr for an isothermal sphere. The ratio of the mean density 
or pressure to the surface value is less than 4 for any negative-index polytrope. On the other hand, 
the density and pressure contrasts between the center and the edge of a critically stable polytrope 
rise from the values shown in Figure ^ as 7 becomes greater than 7 p . In the region denoted "SPS" 
in Figure ||, locally adiabatic, singular polytropic spheres are formally stable against collapse; the 
region in which spheres with large density contrasts are stable against finite perturbations may be 



more limited, as discussed in §5.4. Non-isentropic polytropes may thus provide an explanation for 
the very large density contrasts observed between cores within molecular clouds and the edges of 
the clouds (see Fig. [R] and Curry & McKee 1999). 
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APPENDIX 



A. Standard Nondimensional Variables and Equations 

A set of 'standard' dimensionless variables and associated equations is generally used to 
describe the structure of polytropic gas spheres (e.g., piiandrasekhar 1939 ). Polytropes satisfy 



P(r) = K pP (r)^, (Al) 

where K p is independent of position. The standard equations for isothermal spheres (with 7 p = 1) 
are different from those for other polytropic spheres. An advantage of the variables we use (or 
of the homology variables u and [n + l]v) is that a single set of equations applies to both cases. 
Furthermore, our formulation applies to multiple pressure components as well, whereas the standard 
formulation is restricted to single pressure components. Here we compare the standard isothermal 
and polytropic variables and equations to those introduced in this work. Some of these relations 
have been given previously by Stahler (1983), who focused on the case n > 0; his M n and R n are 
proportional to our \i and A, respectively. 

In the standard equations, the polytropic constant n describes the structure. The relationship 
between n and the polytropic index 7 P is 7 P = (n + l)/n. Isothermal spheres have |n| — > oo; the 
case n = — 1 corresponds to an isobaric cloud. 



A.l. Dimensionless Variables and Structure Equations 

The standard dependent variables are: 

\E' = In — (isothermal spheres) , 
P 



l/n 



(other polytropes), 



(A2) 



(A3) 



where p c is the central density. In this work we introduce a nondimensional mass p in equation 
(|i"o| ) and a nondimensional radius A in equation (|Tl| ) as dependent variables. 



The standard independent variable is a dimensionless radius: 



1 c l 

4vr Gp c 



1-1/2 



(isothermal spheres), 



|n+l| Cc 
4tt Gp c 



-1/2 



(other polytropes); 



(A4) 
(A5) 
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the independent variable in our analysis is x = ^ n (Pc/p)- 

The differential structure equations are derived using the equation of hydrostatic equilibrium 
(|7|). In the standard nondimensional formulation, the equilibrium condition is reflected in the 
Lane-Emden equation, 



-7T — ( i 2 —rr ) = e * (isothermal spheres) (A6) 



1 d f ^ 2 d0\ W + 1| 



e « \ e W = "Vrr r (other ^ 

The structure equations used in this work are equations (^), (PH), and, for multi-pressure 
polytropes, equation (p^). 

The boundary conditions for integration are ^ = 0, d^Sf/dt; = at £ = for isothermal spheres, 
and 9 = 1, d9 /d£ = at £ = for other polytropes. In our formulation, /i = and A = at x = 0. 
Note that both £ and \ vanish at the origin. 



A. 2. Dimensional Quantities 



Equation (A3) states p oc 9 n ; together with the ideal gas law, this proportionality implies 



that the temperature T oc 9, so that c 



Chandrasekhar (1939) and Viala and Horedt 



(1974) identify the following additional proportionalities for polytropic structures: ^ oc distance, 
9' oc temperature gradient, 9 n+l oc pressure, and (— \n + l\/[n + 1])£ 2 #' oc M(£). 

In the standard problem, the dimensional radius can be determined using equations (A4) and 



(A5) together with the equation of hydrostatic equilibrium; alternate forms are: 



R 



GM 



R 



We may write the pressure as: 



GM 



n + 1 cf. 



d9 



R 



GM X 



c 2 fj, 



(isothermal spheres), 

-l 

(other polytropes), 
(this work). 



(A8) 

(A9) 
(A10) 

(All) 
(A12) 

M 2Q3^ ( thlS WOTk J- ( Al3 ) 

We have chosen to express R and P in terms of M{R) and c(R) for our variables, since in our 
approach the dependent variables depend only on properties at the surface. Alternatively, one can 



P 



1 



4vrM 2 G 3 



d 4)\-* (isothermal spheres) 



P 



In + 1| 



4tt M 2 G 3 



inn+l 



P 



dU 

(this work). 



(other polytropes), 
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express the results in our variables in terms of the central sound speed; for a single component 
polytrope, the required relation is c 2 = c 2 exp[— (7 P — 1)%]. 



Finally, the mass can be expressed as: 
M 



1 c 3 ^ 2 cM 



M 



(4vr)V2 G 3/2 p l/^ ^ (isothermal spheres), (A14) 
l« + l| 5/2 c 3 . 2 d6 , s IK s 

(47r) i /2(n+1)G 3 /2/) i/^ -it (other p ol y tro p-)> ( A15 ) 



c 3 



M = G z/ 2p l/2 V ( tMs WOrk )- ( A16 ) 

If desired, the last result can be expressed in terms of the central density by using the relation 
p = p c exp(-x). 



A. 3. Conversion Among Dependent Variables 

Isothermal Case: 



Since ^ = %■> we find 



and 



f = (4vr) 1 / 2 Ae x/2 , (A19) 
e-*/ 2 . (A20) 



1 /i _ x/2 



^ (4vr) 1 /2 A 2 



Polytropic Case: 

Since 6* = exp(— x/n), we find 



(4vr) 1 /2 
and 

(4tt)V2 



~~ (4vr)V2 

A = ^±^ ^ (n " 1)/2 , (A22) 



n + 1 \ L ' Z 



d9 In+ll 1 / 2 



(47r)V2( n + i) A 2 



ii p -7pX/2 



(A24) 
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B. DERIVATION OF EQUATIONS FOR THE PERTURBED STRUCTURE 

B.l. Perturbations in p and A 

To assess the stability of multi-pressure polytropes, we must determine the structure they 
assume after a perturbation. We consider Lagrangian perturbations so that SM = 0. For each 
pressure component we have Pi = K p ip^ 1 . Prior to the perturbation, K p i is constant in space; 
after the perturbation, K p i will generally have a different value, and it may depend on position. In 
general, we find 

5lnPi = 6hiK pi -7 pi 6x, (Bl) 

where we have used the fact that p c is defined to be the central density of the unperturbed sphere, 
so that 5\np = —5\- The perturbation in the total pressure is then 



Pi 



5 In K 



pi 



(B2) 



We can evaluate 5x from the fact that, at constant mass, the definition of p implies that 
p oc p 2 P~ 3 / 2 . Varying this relation and combining it with equation (|B2|), we find 



and 



$X 



5lnP 



4-3 7p 
2 

4-3 7p 



Pi 



lp Sln t i + 2j2[^)s\nK ] 



(B3) 



(B4) 



We now apply Lagrangian variations to the structure equations. Starting with equation (16), 
we have 



d&X 
dx 



<$7jp 

7 P 



5 In 



A* 
p? 



(B5) 



It is convenient to insert equation ( Jl6| ) into equation ( p0[ ) before carrying out its variation: 

dlnp = dlnM- -(4-3j p )dx- (B6) 

Varying equations (|B6|) and @), and using equation (p5|) , we obtain the equations for dS In p/dx 
(eq. ^) and d5lii\/dx (W^) i n the text. 



B.2. Evaluation of 5^f p for adiabatic perturbations 

To complete the set of equations for the perturbed structure, we must determine the 
perturbation in the polytropic index, 5^ p . Recall that a perturbed locally adiabatic polytrope 
is not itself a polytrope unless it is isentropic (§2.2). As a result, we must allow for a variation in 
K p i when we evaluate the polytropic index: 

((MP ^(Pi\( din Kpj ^ 



-37- 



from equation (Bib. (In equilibrium, the cloud is a polytrope, K p i is independent of position, and 



7 P reduces to J2 lpi(Pi/P) as it must — see eq. [fl5| Varying equation (B7), we obtain 



$lv 



dpi - i P )s in K pi - ipidpi -i P )$x 



(B8) 



with the aid of equation (|B1| ). We have used the fact that K p i is constant (d In K p i/dx = 0) in the 
initial equilibrium state in order to simplify this result. 

For a locally adiabatic pressure component £, we have Pg = Kgp^ 1 , where the entropy parameter 
Ki is constant during a perturbation. As a result, we have 5\aPg = 7^<51np = —jg5x, so that 



5lnK pe = (7 p£ -je)Sx, 



from equation (|bT|). As a result, such a component contributes a term 



(it ~ l P e) 



dSx 
dx 



liilpe ~ l P )Sx 



(B9) 



(BIO) 



to Sj p . 



For a globally adiabatic pressure component, it is the total entropy of the cloud that remains 
constant during a perturbation. In order to be distinct from the locally adiabatic case, the 
component must be non-isentropic (7^ 7^ 7po)- Redistribution of energy in the cloud during the 
perturbation maintains the polytropic form, P g oc p^ 9 , but changes the value of K pg (the entropy 
parameter K g changes as well). The change in K vg needed to maintain constant entropy will be 
discussed in §6 below. All that we need to note here is that since K pg is independent of position 
for a globally adiabatic pressure component, such a component contributes a term 



(5j 



y(l P g - lp)(S~LnK pg - lpg 5 X ) 



(Bll) 



to (*>7p from equation (p8[). 



In a multi-pressure polytrope, both types of pressure components may be present. For 
simplicity, we shall assume that only one globally adiabatic component is present; the generalization 
to more than one such component is straightforward. We find it convenient to introduce a new 
adiabatic index 



E 



n + 



p 



7p 9 > 



(B12) 



which is the weighted mean of the adiabatic indexes for the locally adiabatic components I and the 
polytropic index for the globally adiabatic component g. This is equivalent to 



7 P + E 



(7i-7 P i)- 



A related quantity that we shall need is 

r ' = E ( jr) 7i(7 P " Ipi) + y-^J Ipgilp ~ Ipg)- 



(B13) 



(B14) 



iei 
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In terms of these quantities, we find 



4-3r 



(B15) 



from equations ( |B3| ) and ( |B4| ). 

Evaluating 5lnj p from equation (|B8|), we then obtain equation (Q) in the text. 



ENTROPY FOR GLOBALLY ADIABATIC PRESSURE COMPONENTS 



C.l. Determination of the Entropy 

Our objective is to evaluate the entropy parameter -Kg, re f (eq. [Q|) for a globally adiabatic 
pressure component g for an arbitrary non-uniform cloud. To do this, we generalize McKee & 
Zweibel's (1995) treatment of Alfven waves to arbitrary values of 7 g and 7 P9 . The equation for the 
internal energy of an ideal gas is 



dUg 

dt 



+ V • (n 9 v + q ff ) + P g V ■ v = «5, 



9' 



(CI) 



where q 5 is the heat flux and S g is the source term for the globally adiabatic pressure component 
(e.g., McKee et al 1987). (Note that it is possible for S g to vanish in our model even in the presence 
of heating and cooling, provided the effects of the heating and cooling are represented by internal 
degrees of freedom.) The total internal energy in the cloud associated with globally adiabatic 
pressure component is U g = J u g dV, where the integral extends over the entire cloud. The rate of 
change of this energy is 



dUg 

dt 



dt 



dV + 



u g v ■ dS. 



Eliminating du g /dt with equation (|C1|) and using the divergence theorem, this becomes 



dUg 
dt 



Jl-P g X7 -v + S g ]dV - J q 3 -dS, 



(C2) 



(C3) 



where dS is an element of the surface bounding the system. For an adiabatic process, the last two 
terms vanish. Note that we do not assume that the heat flux is zero inside the cloud, only at the 
surface. We do assume that there is no exchange of heat with the other pressure components; for 
example, we do not allow for the possibility that energy would be transferred from the waves to the 
thermal motions of the gas by wave damping. The equation of continuity then allows us to rewrite 
equation ( |C3| ) as 



dt 



(C4) 



where we have included the argument in K vg {t) to emphasize that it may change with time. Since 
the internal energy is given by 



U a 



1 



19 



P g dV 



1g 



^- / P ^- l dM, 



(C5) 
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we can eliminate the unknown K pg from equation flC4|), 

The internal energy can also be expressed in terms of the mean square sound speed ((? g ), 

u '-^=rJ{%) ilts ^i M ®- <C7) 



As a result, equation (C6) can be re-expressed as 



dt \-,„-l)dt \K„) 
Integrating this and evaluating the integration constant in the reference state, we find 

«3Sf _1 = K^{c? g )-^-^\ (C9) 

which is the desired result. 

Next, we express the entropy parameter for a globally adiabatic component in a form similar 
to that for a locally adiabatic component, {Kg). We define the mean entropy parameter as 

( 7p g -i A 1 , f 2Ea_zL\ 

(K 9 )\^-^J =— K} 1P9 ~ l9) dM. (CIO) 

Note that in order to keep the notation simple, we have deviated from our convention of using (x) 
to mean M J x dM; we made a similar exception for (Kg). From dimensional analysis (see eq. 
|C^ ) or by direct evaluation, we obtain 

K^~ l = K^ g - l c- 2 ^-^. (Cll) 

Inserting this into equation ( CIO ), we find that (K g ) is just the entropy in the reference state, 

(K g ) = K gM . (C12) 



Equation ( |C9|) relates the entropy parameter K gtTe { = (K g ) to the polytropic parameter K pg 
in terms of (c 2 ), which is proportional to the mean temperature of the pressure component. In 
the particular case of an isentropic component, we have simply (K g ) = K pg , as was clear from 
the definitions. For the globally adiabatic case of greatest interest to us, Alfven waves, we have 
7 S = 3/2 and 7 pg = 1/2, so that (K g ) = (c 2 g ) 2 /K pg . With the aid of equation QC9D , we can express 
the entropy parameter in terms of directly observable quantities if we choose the pressure in the 
reference state to be the surface pressure P gs : 



(K g ) = -f ^ , (C13) 

( T9-TP9 ] 

where ^ g = (c 2 g )/c 2 gs is usually of order unity, and we have included the subscript "s" to emphasize 
that the quantity is to be evaluated at the surface of the cloud. 
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C.2. Evaluation of 5lnK pg at Constant Entropy 



To evaluate (57 P ) 9 in equation (Bll), we need to know how K pg changes as the result of a 
globally adiabatic perturbation. Variation of equation ( |C9|) at constant entropy gives 

(C14) 



(j g - 1)5 In K pg = (j g - -f pg )Sln(c 2 ). 



Now, with the aid of equation (Bl) we find 



Sc 2 Q = ^(5lnPg-5\np) = cl [SlnK pg - ( lpg - . 



Noting that b\nK p „ is constant inside the cloud, we obtain 



5\n(c 2 ,) =5hiK ] 



7 Pg ~ 1 f 9 



Inserting this result into equation ( C14 ), we then find 

2(7,-7 P9 )/(^)c^M 



5 In K 



P9 



(C15) 



(C16) 



(C17) 



with the aid of equation ( B15| ) for <5£. 
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Fig. 1. — Structure of a spatially isothermal polytrope (7 P = 1) in the fi — A plane. The critical 
density ratios for an isentropic sphere (7 = 1) and a globally adiabatic sphere with 7 = 5/3 are 
indicated. 

Fig. 2. — Structure of a spatially isothermal polytrope (7 P = 1) in the P s — V plane. The critical 
point for 7 = 1 occurs at the maximum value of the surface pressure P s . 

Fig. 3. — The central-to-surface density and pressure ratios for isentropic polytropes (7 = j p ). 
The numerical value is for the isothermal case. 

Fig. 4. — The mean-to-surface density and pressure ratios for isentropic polytropes (7 = 7 P ). The 
numerical value is for the isothermal case. 

Fig. 5. — Characteristics of locally adiabatic polytropes in the 7 — 7 P plane. We do not consider 
polytropes in the region 7 < 7 P since they are convectively unstable. Above 7 = 4/3, the polytropes 
are always stable, even for arbitrarily large values of x = ln(/? c //?s)- Stellar polytropes have 7 > 4/3 
so that they are stable, and in addition have 7 P > 6/5 so that they can have P s = at finite radius; 
as a result, they can be used as models for stable stars. Stable, finite-sized polytropes in the rest 
of the parameter space with 7 > 7 P must be confined by the pressure of an ambient medium. 
The locally adiabatic pressure components in interstellar clouds have 7, < 4/3. For polytropes in 
the region labelled "Finite Xa", the polytropes become unstable for a central-to-surface density 
contrast greater than exp(x CT )- In the region labelled "SPS" (singular polytropic sphere), polytropes 
are formally stable even for an infinite density contrast (x = 00). These regions are separated by 
the curve labelled 700, which is given in equation fl69|). It is possible that clouds in part of the SPS 
region, particularly near the 700 curve at small 7 P , are unstable to finite perturbations (§ 5.4. 3] ). 



Fig. 6. — The value of p CT , which determines the critical mass, is plotted for isentropic polytropes 
as a function of 7 P . Also shown is the value of [i for singular polytropic spheres (labeled SPS). For 
polytropes with a critical point, we have // cr (isentropic) > fi cr ^ ^(SPS). 

Fig. 7. — The same as Figure [f| except that the dimensionless critical mass of the isentropic 
polytrope and the dimensionless mass of the singular polytropic sphere are given in terms of the 
rms value of c 2 instead of the surface value. 

Fig. 8. — The central-to-surface pressure ratio for globally adiabatic polytropes at the critical 
point. Since such polytropes cannot be convectively unstable, it is possible to have stable clouds 
with 7 < 7 P . Globally adiabatic polytropes with 7 > 4/3 are unstable to core collapse for 7 p < 6/5. 
The point at 7 P = 0.5 and 7 = 3/2 is appropriate for small amplitude Alfven waves. 

Fig. 9. — The same as Figure ||, except that the curves represent the central-to-surface density 
ratio at the critical point. 

Fig. 10. — The density distribution of a multi-pressure polytrope that has Pb — P w — 10Pth> 
similar to the values observed by Crutcher (1999) in a sample of molecular clouds. The density 
is normalized to the density at the cloud surface, p s , and the radius is normalized to the value 
of (Cg/GPs) 1 / 2 = c s /(Gp s ) 1 / 2 at the cloud surface; thus, at the surface, the abscissa = A. We 
have assumed that 7 Pj # = 1.0 and that the cloud is close to, but not at, the critical state. For 
comparison, we have portrayed the density distribution of a Bonnor-Ebert sphere with the same 
surface pressure and density. The multi-pressure polytrope has a much larger density drop than 
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the Bonnor-Ebert sphere because of the stabilizing effects of the non-isentropic magnetic pressure; 
the mean-to-surface density drops for the two cases are about the same, however. 
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